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CHAPTER  I 


BACKGROUND 

INTRODUCTION 

The  problem  of  accurately  planning  the  accomplishment  of 
defined  objectives  is  extremely  complex.  Within  the  military  environ- 
ment, planning  support  activities  to  achieve  specified  combat  readi- 
ness objectives  is  especially  complicated.  Such  planning  must  be 
conducted  subject  to  numerous  constraints.  These  constraints 
restrict  the  alternative  courses  of  action  available  to  planners  and 
frequently  cause  the  cost  of  achieving  stated  objectives  to  be  very 
high. 

In  Israel,  this  problem  is  compounded  considerably  by  the 
imposition  of  two  overriding  constraints.  These  constraints  are 
unique  to  the  Israeli  military  forces  and  seriously  limit  the  range  of 
feasible  alternatives  to  those  in  which  cost  minimization  is  not  a 
prime  consideratim.  The  result  is  that  optimal  solutions  tend  to  be 
very  expensive.  A planning  process  in  which  cost  considerations  are 
of  only  secondary  importance  would  not  be  acceptable  to  nations  that 
have  less  severe  military  constraints. 
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The  first  of  these  overriding  constraints  is  imposed  by  the 
ever  present  threat  of  a military  invasion  of  Israel  by  its  neighboring 
countries.  This  threat  makes  it  imperative  that  Israeli  military 
forces  are  maintained  at  a high  level  of  combat  readiness  at  all  times. 

The  second  major  constraint  is  a result  of  limited  national 
resources.  As  a result  of  these  resource  limitations,  Israel  is 
forced  to  rely  heavily  on  foreign  suppliers  of  war  materiel.  This 
dependency  on  external  sources  of  supply  constitutes  a serious  con- 
straint on  defense  planning.  The  significance  of  this  constraint  is 
heightened  when  one  considers  the  impact  foreign  embargoes  would 
have  on  the  Israeli  supply  system.  The  possibility  of  embargoes  is 
very  real,  as  evidenced  by  the  French  embargo  on  military  supplies 
that  occurred  in  1967,  and  the  United  States'  "reassessment  policy" 
that  was  implemented  by  the  Ford  administration  in  1975. 

This  background  information  is  presented  to  emphasize  the 
fact  that  in  many  instances  military  preparedness  planning  is  largely 
predetermined  by  a well-defined  national  policy  that  must  be  achieved 
subject  to  very  binding  constraints.  Thus,  the  Israeli  military 
planner's  range  of  feasible  alternatives  is  quite  restricted. 

GENERAL  SCOPE  AND  JUSTIFICATION 

Maintenance  activities  play  a major  role  in  the  combat  readi- 
ness posture  of  the  Israeli  Air  Force.  Without  a strong  maintenance 
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function,  the  efiectiveneee  of  the  fighting  forces  would  be  greatly 
impaired.  Within  the  Israeli  Air  Force,  the  maintenance  function  is 
divided  into  two  well-defined  subsystems.  These  two  subsystems 
may  be  characterized  by: 

a.  base  level  maintenance  activities 

b.  depot  level  maintenance  activities. 

The  two  subsystems  are  interrelated  through  the  Air  Force  Central 
Inventory  Management  System  (CIMS). 

The  Central  Inventory  Management  System  is  a computerized 
management  information  system  that  controls  all  the  inventory  activ- 
ities t?»t  take  place  within  the  Air  Force.  Repairable  items  that  can- 
not be  repaired  by  base  level  maintenance  are  sent  immediately  to  the 
central  depot.  The  Central  Inventory  Management  System  determines 
the  distribution  of  serviceable  items  from  the  central  supply  base  to 
the  various  using  organizations. 

It  should  be  emphasized  that  the  1.  A.F.  supply  system  is  a 
"push"  system.  In  a "push"  system,  serviceable  items  are  sent  to 
using  organizations  according  to  criteria  that  are  establisned  by  the 
item  managers  in  response  to  total  provisioning  requirements.  This 
system  differs  considerably  from  a "pull"  system  in  which  service- 
able items  are  sent  to  using  organizations  only  upon  receipt  of  a 
requeet  from  the  user. 
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Since  the  CIMS  controls  both  base  level  and  depot  level 
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I nnaintenance  activities,  it  contains  a data  base  that  provides 

researchers  access  to  a vast  amount  of  information.  This  study 
relies  heavily  on  data  extracted  from  the  central  inventory  manage- 
ment system. 

At  the  depot  level,  maintenance  activities  can  be  further 
broken  down  into  three  separate  categories: 

a.  overhaul  and  repair  of  subsystems  and  items 
• b.  overhaul  and  repair  of  main  systems 

c.  overhaul  and  repair  of  aircraft 

I This  research  focuses  on  those  activities  that  fall  in  the  first  cate- 

^ gory,  overhaul  and  repair  of  subsystems  and  items  in  the  depot. 

STATEMENT  OF  THE  PROBLEM 

One  factor  that  greatly  influences  the  ability  of  the  Israeli 
Air  Force  to  achieve  its  desired  preparedness  posture  is  the  exist- 
ence of  a satisfactory  level  of  serviceable  inventory.  Maintaining  a 
satisfactory  inventory  of  serviceable  items  can  be  accomplished  in 

t 

two  different  ways: 

a.  through  the  procurement  process 

i 

b.  through  the  recovery  of  repairable  items. 

The  dependency  of J4rmel  on  external  sources  of  supply,  as  well  as 
budgetary  cons^rations  that  restrict  the  quantities  that  can  be 
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purchased,  forces  the  Air  Force  to  rely  heavily  on  an  efficient 
recovery  system. 

The  capacity  of  the  central  maintenance  depot  is  restricted 
by  the  limited  availability  of  manpower,  equipment  and  monetary 
resources.  This  limited  capacity  imposes  a major  constraint  on 
planning  for  a recovery  system. 

The  purpose  of  this  research  is  to  develop  a mathematical 
model  that  can  be  embedded  within  an  existing  computerized  manage- 
ment information  system  to  provide  depot  level  maintenance  managers 
with  periodic  reports  that  specify: 

a.  The  quantity  of  each  item  that  must  be  recovered  during 

the  next  period. 

b.  The  priority,  or  relative  essentiality,  of  each  repairable 

item. 

More  specifically,  this  research  is  intended  to  provide  a 
means  whereby  each  depot  department  supervisor  can  be  made  aware 
of  the  quantity  of  each  item  to  be  recovered,  as  well  as  the  priority 
that  should  be  given  each  item.  The  priority  of  repair  reflects  the 
relative  essentiality  of  the  item  when  the  total  needs  of  the  Air  Force 
are  considered. 
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BASIC  ASSUMPTIONS 


1.  The  proposed  model  requires  as  inputs  only  data  that 
are  available  in  the  existing  central  inventory  management  system. 

It  is  assumed  that  the  data  are  arranged  in  such  a manner  that  it  will 
be  possible  to  extract  required  data  elements  in  any  stratified  form. 
An  additional  assumption  is  that  it  is  possible  to  extract  historical 
data  that  are  necessary  to  empirically  estimate  the  various  model 
parameters. 

2.  It  is  not  the  intention  of  this  study  to  interfere  with  the 
internal  work  scheduling  activities  of  each  department.  Work 
schedulmg  will  remain  the  responsibility  of  the  production  and  main- 
tenance control  supervisors  within  each  department.  The  intention  of 
the  proposed  model  is  to  enhance  the  capability  of  each  department 

to  meet  serviceable  item  requirements. 

3.  The  model  is  designed  to  conform  to  the  declared  policy 
of  I.  A.F.  Headquarters.  This  policy  states  that: 

a.  The  Air  Force  must  be  prepared  to  sustain  combat 
operations  for  a period  of  X days.  This  requirement  defines  the 
levels  at  which  the  various  serviceable  items  must  be  maintained. 

b.  The  consumption  of  serviceables  in  X fighting  days 
Is  equivalent  to  the  consumption  in  Y peacetime  days. 

c.  Each  major  system  has  the  same  essentiality  to  the 

Air  Force.  For  example,  no  difference  in  essentiality  is  specified  ‘ 
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among  fighter  aircraft,  transport  aircraft,  and  missiles. 

d.  Because  of  the  high  level  of  readiness  that  the  Air 
Force  must  maintain,  every  effort  must  be  made  to  minimize  the  not 
operationally  ready- supply  (NORS)  rate  as  much  as  possible,  with 
relatively  little  attention  being  given  to  marginal  cost  considerations. 

4.  The  proposed  model  does  not  address  cost  versus  effec- 
tiveness analysis.  Exclusion  of  this  kind  of  analysis  is  due  to  the 
previously  described  Air  Force  policy.  There  is  no  doubt  that  in  the 
absence  of  this  policy  the  development  of  a model  for  determining 
recovery  requirements  would  definitely  consider  the  marginal  cost 
of  item  recovery.  However,  I.  A.F.  policy  precludes  estimation  of 
the  cost  of  a nonoperational  system  due  to  a shortage  of  recoverable 
items  (NORS). 

RESEARCH  OBJECTIVE 

The  objective  of  this  research  is  to  develop  a mathematical 
model  that  will: 

a.  determine  the  quantity  of  each  item  to  be  recovered  dur- 
ing each  period  (a  period  is  defined  as  the  time  between  two  succes- 
sive computer -prepared  reports). 

b.  determine  the  priority  that  should  be  asjigned  each 
recoverable  item  when  scheduling  repair  activities. 


The  deeirability  of  being  able  to  accurately  forecast  future 


requirements  is  obvious. 

Many  studies  of  Air  Force  problems  point  out  one 
salient  fact:  An  improvement  in  demand  prediction 
could  result  in  important  savings,  perhaps  more  impor- 
tant than  could  any  other  single  change;  even  a small 
improvement  is  likely  to  yield  significant  savings  [l:l]. 

Increased  accuracy  in  demand  forecasts  will  not  only  result  in  mone- 
tary savings,  but  will  also  result  in  improved  operational  readiness 
and  effectiveness. 

Predictions  of  future  requirements  can  also  affect  planned 
military  operations.  For  example,  forecasted  requirements  that  are 
inaccurately  high  might  result  in  a plarned  operating  being  reduced  in 
strength  or,  perhaps,  even  cancelled. 

In  determining  the  quantity  of  each  item  to  be  recovered, 
several  factors  may  be  relevant.  These  factors  include: 

a.  estimates  of  future  consumption 

b.  total  quantities  of  serviceable  items 

c.  total  quantities  of  repairable  items 

d.  the  level  of  serviceable  items  necessary  to  meet  combat 

readiness  requirements. 

Determination  of  the  priority  of  each  recoverable  item  is 
accomplished  in  accord  with  1.  A.F.  policy  which  states  that  each 
major  system  is  as  essential  as  any  other.  (The  proposed  model  may 
be  easily  modified  to  reflect  any  changes  in  this  policy. ) In 
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determining  a method  for  eatabliahing  priorities,  several  factors 
must  be  considered.  These  factors  include: 

a.  the  number  of  using  bases  for  each  item 

b.  estimates  of  future  consumption 

c.  available  quantities  of  recoverable  items 

d.  available  quantities  of  serviceable  items 

It  was  concluded  that  in  order  to  achieve  the  main  research 
objective  it  would  be  necessary  first  to  accomplish  the  following: 

a.  Identify  the  factors  that  affect  the  model. 

b.  Find  the  relative  weight  of  each  factor  in  the  model. 

RESEARCH  QUESTION 

The  nnajor  question  addressed  by  this  study  may  be  stated 
as  follows: 

Can  a mathematical  model  be  developed  which  will,  when 
embedded  within  an  existing  computerized  management  information 
system,  provide  department  managers  at  the  Israeli  Air  Force 
central  maintenance  depot  with  periodic  reports  that  specify: 

a.  the  quantity  of  each  repairable  item  that  must  be 

recovered  during  the  next  period,  and 

b.  the  priority,  or  relative  essentiality,  of  each  repair- 
able item. 
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CHAPTER  U 

METHODOLOGY 

GENERAL  OVERVIEW  AND  DESCRIPTION 
OF  THE  MODEL 

A*  stated  previously  the  objective  of  this  research  is  to 
develop  a mathematical  model  that  will:  (a)  determine  the  quantity  of 
each  item  to  be  recovered  during  a certain  period,  and  (b)  determine 
the  priority  that  should  be  assigned  to  each  recoverable  item. 

The  variables  that  are  taken  into  consideration  by  the  model 
can  be  classified  into  two  categories: 

a.  Deterministic  variables 

b.  Stochastic  variables 

Deterministic  Variables --Deterministic  variables  are  those  variables 
whose  values  can  be  obtained  directly  from  the  CIMS.  Values  for 
these  variables  are  entered  into  the  model  as  they  appear  in  the  data 
base  of  the  CIMS.  For  eouunple:  number  of  repairable  items  on  hand. 

Stochastic  Variables--Values  for  the  stochastic  variables  represent- 
ing item  demand  are  predicted  using  a mathematical  forecasting 
Ha  model  that  will  be  discussed  later.  Predicted  values  of  this  variable 
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are  entered  into  the  model  in  addition  to  values  for  the  deterministic 
variables. 

In  our  assessment,  the  ability  of  the  proposed  model  to 
accurately  forecast  item  demand  is  a major  factor  in  effectively 
determining  the  quantity  of  each  item  to  be  recovered  each  period 
as  well  as  the  priority  of  each  item. 

DEFINITION  OF  TERMS 

1.  Depot  Level  Recoverable  ltem--This  is  an  item  that  can 
be  repaired/overhauled  only  at  the  depot  level.  Repair  of  a depot 
level  recoverable  item  is  not  authorized  at  the  base  level.  Such  items 
are  sent  inimediately  to  the  depot.  This  item  is  identified  in  the  com- 
puter files  by  code  No.  81. 

2.  Recoverable  Item  Warehouse  (RIW)--Thi8  is  a warehouse 
that  operates  within  the  jurisdiction  of  the  depot.  The  warehouse  func- 
tions as  a central  point  to  which  all  recoverable  items  are  forwarded 
from  the  operational  units. 

3.  Depot  Shop  (DS)--This  is  a workshop  within  the  depot  that 
is  responsible  for  a particular  class  of  recoverable  items. 

4.  Clas8--A  group  of  items  whose  recovery  requires  simi- 
lar equipment  and  knowledge,  e.g.  wheel  and  brake  components. 
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DEFINITION  OF  VARIABLES 

1.  Consumption  of  item  j (Cj)--The  quantity  of  the  item 
actually  consumed  by  users.  This  variable  is  divided  into  four  cate- 
gories: 

a.  Monthly  Consumption  (C^j)--The  quantity  of  the 
item  consumed  in  month  i (i-1,  2. . . 36).  This  information  is  stored 
in  current  data  files  for  the  last  36  months. 

b.  Yearly  Consumption  (C^j)--The  quantity  of  item  j 
consumed  during  the  year«  where  k is  the  number  of  years  the 
item  has  been  used  in  the  Air  Force.  Yearly  consumption  data  arc 
available  in  historic  files  for  every  year  since  1969. 

c.  Predicted  Consumption  (Cpj)--The  predicted  con- 
sumption of  item  j for  the  next  period. 

d.  Predicted  Quantity  Level  (Qpj)--The  predicted  quant- 
ity of  item  j that  is  consumed  in  Y peacetime  days  (which  is  equivalent 
to  X combat  days). 

2.  Total  Serviceable  Quantity  (Q|j)--The  total  serviceable 
quantity  of  item  j in  the  Air  Force. 

3.  Total  Repairable  Quantity  (Q^j)--The  total  repairable 
quantity  of  item  j in  the  Air  Force. 

4.  Priority  (Pj)--The  essentiality  of  item  j relative  to  other 
items  repaired  in  the  same  shop  during  a certain  period. 

' ' '• 
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5.  Quantity  to  be  Recovered  (Xgj)--The  quantity  of  item  j 
that  muat  be  recovered  during  the  next  period. 

6.  Diaperaion  Ratio  (Dj)--The  ratio  of  the  total  number  of 
aerviceable  itema  to  the  number  of  uaera. 

POPULATION  OF  INTEREST 

Thia  reaearch  focuaea  upon  a population  that  ia  characterized 

by; 

a.  recoverable  itema  with  code  81,  and 

b.  itema  that  have  achieved  a a table  trend  of  demand 
behavior.  The  aize  of  thia  population  ia  eatimated  to  be  approximately 
10,000  itema. 

DATA  REQUIREMENTS 

Aa  previoualy  atated,  the  population  of  intereat  ia  aaaumed 
to  have  achieved  a atable  trend  of  demand  behavior.  Conaequently, 
the  data  pertain  to  thoae  itema  that  have  been  in  operation  for  at  leaat 
four  yeara. 

The  fact  that  the  model  containa  both  determiniatic  and 
atochaatic  variablea  requirea  that  data  be  extracted  from  two  different 
aourcea: 

a.  Valuea  for  the  determiniatic  variablea  were  extracted 
from  current  data  fUea. 
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b.  The  data  used  to  determine  values  for  the  stochastic  vari- 
able (the  predicted  consumption)  were  extracted  from  historical  files. 

The  following  table  defines  the  variables  according  to  their 
units  of  measure  and  measurement  scale. 


Table  2-1  (continued) 


demand  (see  "General  Overview  and  Description  of  the  Model"),  and 

b.  research  culminating  in  the  development  of  the  desired 


model  whose  output  will  be  used  by  the  depot  manager. 

Three  mathematical  forecasting  techniques  were  assessed 
in  an  effort  to  determine  which  was  able  to  most  accurately  predict 
future  item  demand.  The  forecasting  methods  evaluated  were: 

a.  Multiple-linear  regression 

b.  Moving  average 

c.  1st  order  exponential  smoothing 

d.  2nd  order  exponential  smoothing 

e.  3rd  order  exponential  smoothing. 

The  second  part  of  the  research  was  concerned  with  using 
the  demand  predictions,  together  with  values  for  the  deterministic 
variables  that  were  obtained  from  current  data  files,  to  establish 
decision  rules  which  determine  the  quantity  and  priority  of  items  to 
be  recovered.  This  part  of  the  research  consisted  of  two  steps: 

a.  Finding  the  logical  relationship  between  the  variables. 

b.  Using  simulation  techniques  to  test  these  logical  relation- 
ships when  different  relative  weights  are  assigned  to  the  variables. 

A flow  chart  depicting  the  research  approach  is  presented  in  Figure 

2-1. 
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> VAUDATING  THE  MODEL 

Validation  of  the  model  was  conducted  in  parallel  with  the 
research  approach  and  also  consisted  of  two  phases: 

a.  Validation  of  the  forecasting  method  used. 

b.  Validation  of  the  overall  model. 

The  first  phase  was  accomplished  by  using  confidence  interval  cri- 
teria to  compare  the  five  forecasting  methods.  The  second  phase 
was  divided  into  two  validating  steps: 

« 

a.  Internal  Tests--these  tests  involved  checking  the  logic 
and  consistency  of  the  model;  e.g. , a higher  serviceable  quantity 
(Qaj)  ehould  cause  a lower  quantity  to  be  recovered  (Xj). 

b.  External  Tests- -the  model  was  compared  with  the  real 
world  thru  simulation  techniques;  e.g.,  actual  data  were  provided  as 
inputs  to  the  model  and  the  resultant  output  was  tested  against  the 
criterion  of  being  able  to  achieve  a certain  level  of  serviceable  items. 


CHAPTER  m 


DEVELOPMENT  OF  THE  FORECASTING  MODEL 

INTRODUCTION 


The  objective  of  this  chapter  is  to  compare  several  methods 
used  to  forecast  recoverable  items  demand  requirements.  This 
problem  is  important  due  to  the  excessive  cost  incurred  as  a result 
of  inaccurate  forecasts.  Inaccurate  prediction  of  recoverable  items 
demand  results  in  two  types  of  unnecessary  costs.  Associated  with 
predicting  a requirement  much  greater  than  actually  needed  are  the 
costs  of  acquisition  and  of  storage  or  disposal  of  unused  items.  When 
the  predicted  requirement  is  less  than  the  actual  demand,  excessive 
reordering  costs  are  incurred.  More  importantly,  the  level  of  mis- 
sion readiness  will  be  adversely  affected. 

To  make  a prediction  of  demand  for  recoverable  items,  it  is 
necessary  to  know  what  demand  has  been  in  the  past.  The  data  used 
for  making  forecasts  of  recoverable  items  consumption  is  obtained 
from  items  of  the  same  class.  ^ 

^Class  is  defined  as  a group  of  recoverable  items  which  need 
the  same  knowledge  and  equipment  while  repaired. 


The  performances  of  five  commonly  used  forecasting 
methods  will  be  compared.  These  five  methods  are: 

a.  Moving  Average 

b.  First  order  exponential  smoothing 

c.  Second  order  exponential  smoothing 

d.  Third  order  exponential  smoothing 

e.  Multiple-linear  regression. 

TIME  SERIES  MODEL 

A set  of  ordered  observations  of  a quantitative  variable  taken 
at  successive  points  in  time  is  called  a time  series.  Time  in  terms 
of  years,  months,  days,  or  hours  is  simply  a device  that  enables  one 
to  relate  all  phenomena  to  a set  of  connmon.  stable  reference  points. 

A typical  time  series  may  be  thought  of  as  being  made  up  of 
four  types  of  fluctuations. 

1.  A trend  component  which  may  be  defined  as  long  term 
growth  or  decay. 

2.  A seasonal  component  which  may  be  defined  as  a regu- 
larly recurring  periodic  fluctuation. 

3.  A cyclical  component  which  may  be  defined  as  a wave-like 
fluctuation  about  trend.  The  length  and  the  cycle  are  not  constant  as 
in  tha  sastsonal  component  but  may  vary  from  one  cycle  to  the  next. 

4.  An  erratic  or  random  component  which  is  completely 
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If  the  letters  T,  S,  C,  and  E are  used  to  represent  the  four  ^ 

components  of  trend,  seasonal,  cyclical,  and  erratic  fluctuations, 
respectively,  the  time  series  Y,  may  be  written  in  two  basic  models. 

1.  The  multiplicative  model  in  which  Y = TxSxCxE 

2.  The  additive  model  in  which  Y sT+S  + C+  E 

Moving  Average 

The  moving  average  is  the  simplest  forecast  procedure  for 

1 

time  series  analysis  or  forecasting.  This  method  damps  noise  by 
averaging  several  pieces  of  data.  By  updating  the  average  sequen- 
tially by  adding  in  each  new  piece  of  data  while  dropping  off  the 
oldest,  process  changes  can  be  gradually  tracked.  According  to  this 
procedure  the  current  average  is  the  forecast  for  all  future  periods, 

< 

I but  unfortunately  the  moving  average  lags  considerably  more  abrupt 

I changes  in  the  time  process;  e.  g.  an  n-period  moving  average  will  not 

reflect  completely  a step  change  for  n-periods. 

I<et  represent  the  observations  at  time  t where  t = 1,2,... 

N.  The  number  of  the  most  recent  observation  as  used  in  the  single 
moving  average  is  N.  Then  the  actual  average  of  the  N most  recent 
observations,  taken  at  time  t is 

^ N 

Each  of  the  past  N observations  is  weighted  1/N;  all  the 
older  observations  receive  aero  weight.  For  the  next  period,  the 
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data  includea  the  newest  observation,  X 


and  excludes  the  oldest 


.1 

I 


t+1' 


observation 


M 


t+1 


^ ^t+1  Xt  + . . . Xt.i^.^2 

N 


Thus  the  data  base  remains  N periods  behind  the  forecast  periods, 
deleting  the  oldest  data  found  and  adding  the  newest  observation  to 
the  data  base.  This  accounts  for  the  terminology,  moving  average. 


First  Order  Exponential  Snioothing 

Exponential  smoothing  is  so  named  because  the  weights 
assigned  to  demand  observations  are  constructed  by  raising  a fraction 
to  successively  higher  powers  or  exponents.  Exponential  smoothing 
is  also  known  as  exponential  forecasting,  adaptive  forecasting,  adap- 
tive smoothing  or  geometrical  smoothing. 

Exponential  smoothing  is  a special  kind  of  moving  average 
that  does  not  require  keeping  a long  historical  record  in  the  active 
file  and  thus  reduces  the  data-processing  time  required.  Like  other 
moving  averages,  it  has  a response  to  change  which  does  not  fluctuate 
rapidly,  but  the  rate  of  response  can  be  adjusted  readily.  (This  rate 
of  response  is  explained  in  greater  detail  below. ) Exponential  smooth- 
ing is  deemed  appropriate  if  the  data  follow  the  constant  time  series 
model  as  represented  by  the  equation 

a>  a + 
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where 


ie  the  observation  at  time  t 
a is  a constant 

is  the  random  fluctuation. 

Building  the  Model.  The  fundamental  idea  of  exponential  smoothing  is 
that  any  new  estinnate  of  demand  is  based  on  an  old  estimate  corrected 
for  current  information.  The  estimating  equation  is: 

New  estimate  = (old  estimate) 

+ a (new  demand  - old  estimate)  (1) 

or 

New  estimate  = a (new  demand) 

+ (1-a)  (old  estimate)  (2) 

The  new  estimate  is  a smoothed  value  and  is  equal  to  the  previous 
smoothed  value  plus  a fraction  £ of  the  difference  between  the  new 
observation  and  the  previous  smoothed  value.  Letting  St(x)  equal  the 
new  estimate  or  smoother  value  at  time  ^ a equal  the  smoothing 
constant,  and  equal  the  latest  observation  at  time  ^ equation  (2) 
becomes 

St(x)  « aXt  + (l.ar)St.i(x)  (3) 

The  choice  of  the  smoothing  constant  a is  arbitrary  and  is  covered  in 
depth  later.  The  data-processing  simplicity  is  obvious,  since  only 
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one  number  S(_^(x)  has  to  be  recorded  instead  of  the  actual  demand 
in  each  of  the  past  ^ months.  The  data  for  the  past  N months  are 
combined  to  yield  S^_j(x)  in  the  following  manner: 

St(x)  = aXt  + (l-a)St.i(x) 

= oXt;  + (l-a)[aXt.i  + (l-o)Sj_2(x)] 

= oXt  +flr(l-a)Xt.i  + (l-a)^[aXt.2  + ( 1 - a)  S^_3(x)] 

= oXj  +a!(l-a)Xt.i  + o(l.a)2Xt.2  + ... 

+ a(l-a)"Xt.^  + ...  + (l-o)tXo  (4) 

The  function  S^(x)  is  a linear  combination  of  all  past  observations. 
Since  the  expectation  of  the  function  is  equal  to  the  expectation  of  the 
data,  the  function  can  be  called  an  average.  Thus,  there  is  a new 
way  of  estimating  the  value  of  the  coefficient  in  a constant  model. 

= St(*)  (5) 

Selecting  the  Smoothing  Constant.  The  analysis  of  the  factors  that 
influence  the  choice  of  the  smoothing  constant  (or  the  number  of  the 
periods  of  data  that  are  average)  is  basic  to  the  concept  of  making  a 
satisfactory  compromise  between  stability  of  the  forecasts  and  the 
speed  of  response  to  changes  in  the  model  that  represents  the  data. 

The  value  chosen  for  the  smoothing  constant  a determines 
what  effect  various  data  points  have  upon  the  estimate  of  the  average. 
As  in  the  case  of  the  moving  average,  the  more  past  data  included  in 


the  average,  the  enialler  the  error  in  the  estimate--provided,  of 
course,  that  the  basic  pattern  of  denumd  does  not  change  during  the 
interval.  On  the  other  hand,  if  fewer  past  periods  are  included  in 
the  averaging  process,  the  response  is  faster  to  the  changes  that  do 
occur. 

As  a practical  matter  the  choice  of  or  would  usually  be 
arrived  at  by  trial  and  error  using  a sample  of  actual  past  demand  or 
usage  data.  An  a would  be  tried;  the  errors  could  be  calculated  and 
compared  with  other  a's.  The  value  of  or  minimizing  error  variance 
would  then  be  a good  choice. 

Weighting  past  data.  Recall  that  the  smoothing  constant  controls  the 
number  of  past  observations  that  have  any  effect  on  the  forecast.  The 
smoothing  constant  or  can  take  on  any  value  from  (but  not  including)  0 
to  1.  When  a is  chosen  very  close  to  0,  this  is  essentially  equivalent 
to  using  an  arithmetic  average  over  a long  period  of  time  as  the  best 
estimate  of  the  future  usage  rate.  When  a is  1,  the  formula  reduces 
to  using  the  current  period's  usage  rate  as  a forecast  or  estimator. 
Intermediate  choices  for  or  produce  forecasts  which  give  more  or  less 
emphasis  to  long-run  average  versus  current  usage,  depending  on 
whether  a is  chosen  closer  to  0 or  to  1. 

For  a first  order  exponential  smoothing  model,  the  weight 
attached  to  historical  observations  decreases  exponentially  according 
to  (1-0 )^,  k > 0,  1,  .. .,  n. 
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changes  in  the  time  series  is  a function  of  the  relative  sise  of  a.  The 


lower  a,  the  slower  the  response.  Higher  values  of  a cause  the 
smoothed  average  to  react  quickly- -not  only  to  real  changes  but  also 
to  random  fluctuations. 

Average  age  of  the  data.  One  way  to  define  an  exponential  smoothing 
system  that  is  equivalent  to  an  N-period  moving  average  is  to  say  that 
the  smoothing  constant  is  selected  to  give  the  same  average  age  of  the 
data.  In  a moving  average  scheme,  each  of  the  N most  recent  obser- 
vations is  weighted  equally  by  1/N,  and  all  prior  observations  are 
weighted  zero.  Let  the  average  age  of  the  data  used  in  a moving 
average  be  defined  as 

K - - 2 (o) 

The  average  age  is  the  age  of  each  piece  of  data  used  in  the  average, 
weighted  as  the  data  of  that  age  would  be  weighted.  In  the  exponential 
smoothing  process,  the  weight  given  data  k periods  ago  is  am  » 
where  /3  s 1-a.  Hence,  the  average  of  the  age  data  is 

IT  = Oo  + lo/8  -f  + ... 

• /8/o  . (l-o)/o  (7) 

k»l 
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Therefore,  to  determine  a single  moving  average  (SMA)  or  single 
exponential  smookfaing  (SES)  process  of  equal  average  age  of  data, 
when  a in  SES  of  N in  SMA,  respectively,  is  known,  solve 


1-0 

a 


N-1 

2 


or  a 


2 

N+1 


(8) 


It  is  not  sufficient  to  equate  the  average  ages,  however,  as  the  SMA 
weights  the  data  equally,  and  the  SES  process  places  heavy  weight  on 
the  most  recent  observations  and  less  on  the  older. 
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Initial  Conditions.  Repeating  Equation  (3) 


1 ^ 


it  will  be  noticed  that  exponential  smoothing  always  requires  a prev- 
ious value  of  the  smoothing  function.  When  the  process  is  started, 
there  must  be  some  value  that  can  be  used  as  the  previous  value  S(._^ 
(x).  If  there  are  past  data  at  the  time  one  starts  to  use  exponential 
saooothing,  then  the  best  initial  value  would  be  a simple  average  of  the 
most  recent  N observations:  « M^_^  initially. 

Frequently,  there  is  no  past  data  to  average;  smoothing  starts 
with  the  first  observation.  In  this  case,  a prediction  of  the  average  is 
required.  The  prediction  may  be  what  the  process  is  intended  to  do. 

In  other  cases,  the  prediction  can  be  based  on  similarity  with  other 
processes  that  have  been  observed  for  some  time,  as  in  the  case  of 
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a new  item  added  to  an  inventory. 

Summary  of  Single  Exponential  Smoothing.  Exponential  smoothing  is 
one  of  the  simplest  yet  most  flexible  methods  of  forecasting  demand 
from  past  observations.  It  is  accurate.  It  is  simple,  at  least  by 
comparison  with  any  moving  average  technique,  requiring  only  two 
multiplications  and  one  addition  and  using  only  two  pieces  of  data,  the 
forecast  and  the  observation  for  the  latest  period.  It  is  flexible 
because  the  weight  assigned  to  each  period's  observation  can  be  easily 
changed.  Thus  exponential  smoothing  meets  both  criteria  of  (1)  sim- 
plicity of  computation,  and  (2)  flexibility  to  adjust  the  rate  of 
response. 

The  demand  pattern  must  fit  the  constant  model  represented 
by  equation  (1).  If  the  demand  pattern  exhibits  a linear  trend,  i.  e. , 
the  process  is  changing  at  a steady  rate,  the  forecasts  made  with  SES 
will  lag  behind  the  actual  time  series  observations.  This  lag  can  be 
corrected  by  introducing  a trend  component  into  the  forecasting  equa- 
tion. This  is  accomplished  in  double  exponential  smoothing  (DES). 

Second  Order  Exponential  Smoothing 
The  fundamental  equation  is 

St(x),  2 aaSj(x),  1 + (l-a)Sj^j(x),  2 

30 


Estimates  of  the  coefficients  of  the  first  order  polynomial  (linear) 
model  are  obtained  from  the  first  and  second  order  smoothed  sta- 
tistics. 

S^(z),  1 and  S^(x),  Z,  through  the  following  equations 

fio(t)  = 2 Sj(x),  1 - Sj(x),  2 

bi(t)  *(-^)[St(x),  1 -St(x),  2] 

the  initial  conditions  to  start  smoothing  are: 

Sgix),  1 =X(0)  .(i^)tj(O) 

Sq(x).  2 =X(0)  - 2(i^)  bi(0) 

where  &j(0)  is  an  initial  estimate  of  the  slope  at  t = 0. 

The  general  equation  for  the  forecast  is 

St(T)  = bQ(t)  + bj(t)(T-t) 

Note:  Selecting  the  smoothing  constant  a is  done  by  the  same  proce- 
dure described  in  the  first  order  exponential  smoothing. 

Third  Order  Exponential  Smoothins 
The  fundam«ital  equation  is 

S^:x).  3 - aSj(x),  2 + (l-or)  S^.i(x),  3. 
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Estimates  of  the  coefficient  of  the  second  order  polynomial  (quad- 
ratic) model  are  obtained  from  the  first,  second,  and  third-order 
smoothed  statistics,  S^(x),  1,  S^(x),  2 and  S^(x),  3 through  the  folloDv- 
ing  equations. 


bo(t)  « 3St(x),  l-3S^(x),  2 +St(x),  3 

bi(t)  = S'  [6-5o)  S.(x),l  - (10-8o)S.(x),2  + (4-3a) 

2(l-a)2 

(St{x),3)] 

b,(t)  . SL—  [Sj(x),  1 - 2S.(x),2  + S.(x),  3] 

2(l-or)2 

The  initial  conditions  to  start  the  smoothing  are: 

So(x),  1 = bQ(0)  - (1^)  fei(O)  + ^2(0) 

So(x),2  » 6o(0)  - ii^6i(0)  + 

So{x),3  .^0(0)  ^2(0). 

The  selection  of  the  smoothing  constant-a,  is  similar  to  that 
described  above. 


MULTIPLE -LINEAR  REGRESSION 


Much  of  the  effort  expended  in  scientific  research  is  devoted 
to  the  search  of  a cause  and  effect  relationship  that  may  exist  between 
phenomena.  The  process  of  discovery  in  this  deterministic  (cause 
and  effect)  mode  of  inquiry  involves  learning  by  association.  Events 
are  observed  and  related,  and  on  the  basis  of  this  analysis  a decision 
is  made  whether  the  events  are  related  casually.  The  events 
analyzed  include  a set  of  observed  variations  of  a dependent  variable 
and  a set  of  corresponding  variations  for  each  of  one  or  more  inde- 
pendent variables. 

Regression  and  correlation  are  powerful  statistical  tools  that 
provide  quantitative  expressions  or  models  of  the  moment  or  extent 
to  which  events  are  related  mathematically.  The  application  of  these 
■ tatistical  methods  cannot  offer  proof  of  the  existence  of  a causal 
relationship  between  selected  variables.  Statistical  analysis  of  these 
methods  can,  however,  provide  valuable  information  that  the  analyst 
can  employ  to  support  a judgment  concerning  the  existence  of  a cause 
and  effect  relationship  between  the  variables  selected  for  analysis. 

The  problem  of  making  a prediction  or  an  estimate  of  the 
value  of  a dependent  variable  corresponding  to  a given  value  of  an 
independent  variable  requires  that  a mathematical  expression  in  equa- 
tion form  be  devised  to  express  the  relationship  between  the  relevant 

/ 
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variables.  This  equation  is  an  expression  of  the  functional  relation- 
ship between  the  assumed  dependent  variable  and  one  or  more  inde- 
pendent variables.  The  relationship  in  which  a dependent  variable 
can  be  determined  exactly  from  a selected  value  of  one  or  more 
independent  variables  (XpX2,X3, . . . .x^^)  is  described  by  the  func- 
tional expression 

Y = bjxj  + b2X2  + ...  b^Xjj. 

The  most  appropriate  method  used  to  find  the  coefficients 
(bi,b2. . . • fb^)  of  this  relationship  is  the  least  square  method  (LSM) 
which  minimizes  the  sum  of  the  squared  differences  between  the 

A 

actual  values  of  Y and  its  calculated  value  from  the  functional  expres- 
sion Y.  In  mathematical  expression: 

Min  53  (Y-^)2 
i=l 


The  statistical  significance  of  the  overall  MLR  model  and  the 
coefficients  (b|^,b2t  • . • fb^)  are  derived  using  the  F test  and  the  t test 
respectively  (4:604-633). 


METHODOLOGY 


Source  of  Data 

The  data  used  in  generating  the  forecasts  and  in  comparing 
the  alternative  forecasting  methods  were  actual  usage  data  which 
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were  extracted  from  the  ClMS's  files  (the  computerised  inventory 
system  of  the  LAF). 

Description  of  Data 

As  stated  above,  five  methods  of  forecasting  are  compared 
in  this  study.  These  methods  can  be  divided  into  two  main  categories: 

a.  Multiple -linear  regression  (MLR) 

b.  Time  series  analysis. 

MLR  Data  Requirements --Using  the  MLR  method,  future  demand  is 
forecast  as  a function  of  the  following  variables. 

a.  Level  of  activity  (flying  hours) 

b.  Number  of  major  systems  in  which  the  observed  item  is 
installed. 

c.  The  class  to  which  the  observed  item  belongs. 

A stratified  sample  of  30  items  was  randomly  selected  from  the 
following  three  classes. 


Table  3-1 
Sample  Sices 


Class  of  Components 

Population  Size 
(approximate) 

Engine 

1500 

15 

Wheel  li  Brake 

600 

8 

Communication 

700 

7 
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Definition  of  Variablee 


Yj  - The  observed  demand  of  the  item. 

Xjj  - The  observed  number  of  flying  hours  of  the 

item. 

X2j  - The  number  of  major  systems  in  which  the 

item  is  installed  in. 


^Dlj 


XD2j 


- Is  a dummy  variable  which  defines  the  class  of 
the  item- -will  be  "1"  if  it  belongs  to  the 
"engine"  class  and  "0"  otherwise. 

- Is  a dummy  variable  which  defines  the  class  of 
the  items --will  be  "1"  if  it  belongs  to  "wheel 
& brake"  components  and  "0"  otherwise. 


Table  3-2  illustrates  the  assignment  of  values  to  dummy  variables. 


Table  3-2 

Identification  of  Dummy  Variables 


Dununy 

Class  .J^^iable 

Xdi 

Xd2 

Engine 

1 

0 

Wheel  b Brake 

0 

1 

Communication 

0 

0 

Note:  The  number  of  dummy  variables  is  the  number 
of  classes  minus  one. 


This  model  considers  only  the  independent  variables  which 
were  mentioned  above;  however,  in  our  assessment,  other  indepena- 
ent  variables  such  as  number  of  landings,  type  of  mission,  number 
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of  missions,  may  possibly  influence  the  dependent  variable.  These 
variables  were  not  considered  because  they  were  not  available  from 
our  data  source. 

Time  Series  Data  Requirements --In  the  time  series  analysis,  24 
values  of  observed  demand  were  used  to  forecast  the  future  demand, 
(24  months  of  past  consumption).  These  values  were  obtained  for  one 
item  of  each  class  that  was  mentioned  previously. 

MLR  Analysis  Approach 

The  analysis  of  the  MLR  model  employed  an  existing  library 
program,  MLREG,  which  uses  LSM.  This  program  enabled  us  to: 

1.  Determine  the  coefficients  of  the  following  regression 
equation 

y = '^0‘*’^1*1‘’’^2*2'*'*^D1*D1+^D2*D2 

2.  Test  the  significance  of  overall  model. 

3.  Test  the  significance  of  the  regression  equation  coef- 
ficients. 

4.  Find  the  net  or  marginal  contribution  of  each  variable  to 
the  model's  explanatory  power. 

Time  Series  Analysis  Approach 

The  analysis  of  the  time  series  method  required  the  use  of 
two  conq^teriaed  programs: 
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1.  SLSC/Moving,  R - -for  the  moving  average. 

2.  SI.  LIB/TCAST--£or  the  exponential  smoothing  model. 

The  first  program  used  24  data  points  to  forecast  the  demand 

for  the  25^^  month  using  a moving  average  model  with  N = 2, 4,6,9,  12. 
For  each  of  the  five  moving  averages,  the  mean  absolute  deviation 
(MAD  was  recorded).  The  optimal  moving  average  model  was  the  one 
with  the  minimum  MAD. 

The  second  program  (TCAST)  forecasted  the  demand  in  the 
25^  month  (lead  time  = 1)  with  different  values  for  a (a  = (0.04,  0.08, 
0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.7,  0.8,  0.  9)  and  for  each  of  the 
three  exponential  smoothing  models.  The  preferred  model  (optimum 
a,  type  of  smoothing,  optimal  cyclic  behavior)  was  obtained  by  observ- 
ing the  minimum  MAD. 

MLR  Results  Analysis 

a.  The  significance  of  the  overall  model  was  tested  thru  the 

F statistic  which  was  computed  as  the  ratio  of  the  mean  square  regres- 
sion (MSR)  to  the  mean  square  error  (MSE).  A higher  F statistic 
indicates  a higher  significance  level  of  the  overall  model. 

b.  The  significance  of  each  coefficient  of  the  regression 
equation  was  tested  by  means  of  a t-statistic  which  was  computed  as 
the  ratio  of  the  coefficient  value  to  its  standard  deviation.  A higher 
t-statistic  results  in  a higher  significance  level  of  this  coefficient 
existing  in  the  model. 
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c.  The  localized  or  stepwise  contribution  of  each  variable 
to  the  explanatory  power  of  the  model  is  available  from  the  computer 
output  because  the  program  selects  the  variable  to  be  introduced  next 
on  the  basis  of  its  explanatory  power. 

d.  The  validity  of  the  model  was  tested  by  observing  the 
number  of  actual  values  which  fell  within  the  forecast  confidence 
interval  for  each  forecasted  value. 

The  confidence  interval  for  a single  future  observation  of  the 
process  at  x®  (simple  linear  regression  model)  was  computed  using 
the  following  formula  (for  the  (1-a)  100%  confidence  level). 


C.I. 


A 

= Y 


x0~ 


df  = n-2 


(x0-3E)^ 


where 

Y I Q is  the  forecast  value  at  the  point  x'^ 

^a/2  n-2  appropriate  value  extracted  for  the  t distri- 

' bution  with  a significance  level  of  a and  n-2 

degrees  of  freedom 

dy.j^  is  the  estimate  of  the  standard  deviation  of  dependent 
variable  Y 

n is  the  number  of  observations 

X is  the  mean  of  the  independent  variable. 

This  approach  can  be  expanded  to  more  than  two  dimensions.  When 
the  confidence  interval  is  defined  in  more  dimensions  than  can  be 
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f visualized,  calculations  beccme  easier  when  matrix  algebra  is  used. 

The  interested  reader  can  consult  Applied  Linear  Statistical 
Models  by  Neter  and  'Wasserman,  p.  245,  for  a complete  explanation 
of  the  matrix  approach  to  solving  multiple  linear  regression  models. 

The  confidence  intervals  were  computed  using  an  existing 
computer  program- -MREG,  R.  This  program  splits  the  data  base 
into  two  parts  (according  to  a pre-selected  percentage),  and  builds 
an  MLR  model  from  one  part  euid  computes  confidence  intervals  for 
remaining  data  points. 

Time  Series  Results  Analysis 

The  results  of  the  time  series  forecasting  method  were 
analyzed  for  one  randomly  selected  item  from  each  class  mentioned 
above  as  follows: 

a.  For  the  moving  average  method  the  optimal  N (number  of 
observations  used  to  compute  the  moving  average)  was  extracted  from 
the  computer  output  for  each  item.  The  optimal  N was  the  one  with 
the  minimum  observed  MAD. 

b.  For  the  exponential  smoothing  method  the  preferred  type 
of  smoothing  and  optimal  a were  observed  for  each  item.  These 
values  were  provided  as  an  output  by  the  program  and  in  fact  were 
chosen  according  to  the  minimum  MAD  method. 

c.  The  optimal  parameters  for  these  two  methods  were  then 
compared  as  follows: 
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(1)  The  minimum  MAD  for  each  item  pointed  out  the 
minimal  abaolute  deviation  forecasting  method.  However,  it  should 
be  noted  that  this  procedure  does  not  necessarily  lead  to  the  best 
method,  because  of  the  possibility  of  it  being  inefficient. 

(Z)  The  predictive  efficiency  of  each  method  was  tested 
by  observing  the  number  of  actual  values  which  "fell"  within  the 
confidence  interval  for  each  forecasted  value.  The  model  with  the 
highest  efficiency  level  was  accepted. 

Comparison  Analysis  of  MLR  and  Time  Series 

The  comparison  between  the  MLR  model  and  the  Time  Series 
models  was  accomplished  through  two  methods. 

a.  Con^>ariBon  between  the  MAD  for  each  class  (in  the  time 
series  model  this  class  is  represented  by  one  item). 

b.  Comparison  of  efficiencies  between  the  ''most  efficient" 
time  series  model  and  the  MLR  model.  The  flow  chart  in  Figure  3-1 
describes  the  process  of  analyzing  MLR  and  time  series  results. 

PRESENTATION  AND  INTERPRETATION 
OF  RESULTS 

MLR  Results  Presentation  and  Interpretation 

The  regression  equation  that  was  obtained  from  the  computer 
output  was 
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9 = 3.74811  + 0.  0006lxj  + 0.01223*2 
+ 3.  18313*^^  + 12.  15841  Xjj2 

where 


A 

y 

is  the  forecasted  demand 

*1 

is  the  number  of  flying  hours 

*2 

is  the  number  of  major  systems 

*D1 

is  the  dummy  variable  which  defines 

"engine"  class 

*D2 

is  the  dummy  variable  which  defines 
class. 

"wheel  it  brake" 

1.  SiHnificance  of  the  overall  model- -Table  3-3  illustrates  the  vari- 
ance analysis  of  the  model. 


Table  3-3 

Analysis  of  Variance 


Variation 

D.F. 

Mean 

Square 

F 

Statistics 

Explained  Variation-EV 

803.277496 

4 

200.819374 

6.  802002 

Unexplained  Variation- 
UV 

738.089772 

25 

29.523567 

ToUl  VarUtion-TV 

1541.366669 

29 

Coefficient  of  Determination  R^  = 

0.5211463 

The  critical  value  of  F in  a confidence  level  of  90%  is  2. 18.  The 
appropriate  H.T.  involving  (the  significance  of  the  model)  is 


Fg  = 6.  802002 


r 


■:] 


Hq!R^  = 0 
Hj:R2^0  Fc=2.  18 

4,25,  a = 0.  10 

Fs>  Fq  - reject  Hq  and  conclude  that  the  overall  model  is 
significant  at  a 90%  C.  L. 

2.  Significance  of  the  Coefficients --Table  3-4  presents  the  t statis- 
tics for  the  different  coefficient. 


Table  3-4 
t Statistic 


Var 

Value 

t Statistics 

®0 

3.74811 

0. 64460 

1.708 

Bl 

0. 00061 

0. 06302 

1.708 

B2 

0.01223 

0. 52482 

1.708 

®D1 

3. 18813 

1.05767 

1.708 

»D2 

12.  15841 

3. 68490* 

1.708 

The  appropriate  H.  T.  involving  each  is 


= t 0 

Rejection  of  the  null  hypothesis  will  mean  that  the  marginal  contribu- 
tion of  the  Xj  is  significant  at  a C.  L.  of  90%. 

‘c  * ‘25,0. 10  "1*708 
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The  conclusion  of  this  teat  is  that  only  the  dummy  variable  Xd2 
('wheel  and  brake  class)  is  significant. 

3.  Contribution  of  each  variable  to  the  explained  variation  (E.  V.) 
of  the  model- -Table  3-5  presents  the  contribution  of  each  variable 
to  the  E.V.  of  the  model. 


Table  3-5 

Explained  Variation  for  Regression  Model 


Model 

E.V. 

AE.V. 

AR^ 

y = Bq  + BjXi 

37.528168 

37.528168 

0.0243473 

0. 0243473 

y = Bq'*’®!*!  +®2*2 

398. 889069 

361.360601 

0.2587892 

0. 2344419 

y » BQ  + biXi+B2X2 
®D1*D1  ®D2*D2 

803.  277496 

404.388427 

0.5211463 

0.2623571 

The  conclusions  from  Table  3-5  are: 

a.  The  variable  (flying  hours)  has  a meaningless  contri- 
bution to  the  explanatory  p<mer  of  the  model. 

b.  The  variables  that  contribute  most  are  the  dummy  vari- 
ables which  represent  the  association  of  the  items  to  different  classes. 

c.  The  variable  X2  (number  of  major  systems)  has  a mean- 
ingful contribution  to  the  explanatory  power--in  contrast  to  the  activity 
level  (X|--flying  hours). 

The  conclusion  that  flying  hours  contribute  a meaningless 
amount  to  the  explanatory  power  of  the  model'  is  not  surprising  since 
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several  studies  have  indicated  that  little  or  no  correlation  exists 
between  demand  and  the  activity  level.  For  example,  Denicoff  and 
Haber,  in  a study  conducted  in  19^2,  found  that  there  was  a measur- 
able relationship  between  demand  and  flying  hours  in  only  3%  of  the 
items  observed  (2:9)«  In  a similar  study  conducted  by  Denicoff  it  was 
concluded  that  for  most  items  demand  projections  based  on  flying 
hours  would  be  futile  (3:2). 

Studies  by  the  Rand  Corporation  indicate  that  a major  prob- 
lem in  determining  the  existence  of  a relationship  between  demand 
and  program  activity,  such  as  flying  hours,  may  be  caused  by  a time 
factor.  Demand  for  a particular  item  may,  in  fact,  be  highly  corre- 
lated with  flying  hours  but  if  the  service  life  of  the  part  is  quite  long, 
the  demand  figures  will  not  follow  the  month-to-month  fluctuations  in 
flying  hours  (1:86). 

In  our  assessment,  the  fact  that  the  number  of  major  systems 
was  found  to  be  a contributory  variable  to  the  model  must  be  investi- 
gated more  deeply  in  order  to  be  able  to  draw  correct  conclusions 
about  its  explanatory  power.  The  reason  for  this  assessment  is  that 
no  statistically  significant  conclusions  can  be  drawn  from  the  sample 
that  we  observed,  because  of  its  small  size  in  comparison  to  the 
population  size.  (The  availability  of  data  sources  restricted  the 
sample  size  to  30  observations. ) 

The  fact  that  the  classification  of  the  population  (dummy  vari- 
ables) has  the  most  net  or  marginal  explanatory  power  can  be 
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explained  by  the  fact  that  there  ie  a correlation  between  the  demand 
and  a certain  class  of  items.  However,  we  cannot  conclude  that  a 
statistically  significant  relationship  exists  because  of  the  limitations 
on  sannple  sise. 

4,  MeasurinE  predictive  efficiency  of  the  model- -As  previously 
stated,  the  measurement  of  predictive  efficiency  was  carried  out 
using  the  MKEG  program.  The  data  were  randomly  split  to  80/20 
percent  (80%  used  to  build  the  model).  Table  3-6  presents  the  results. 


Table  3-6 

Confidence  Intervals  for  Regression  Model 


Y 

C.l. 

4 

7.43 

13.29 

10 

8.74 

14.55 

14 

7.25 

15.30 

7 

5.  17 

13.05 

12 

7.33 

13.92 

8 

7.05 

13.  14 

The  large  values  of  the  C.l.  obtained  (approximately  2 times  larger 
than  the  nominal  values)  "ensure"  that  each  predicted  value  "falls"  in 
the  interval.  Therefore  we  can  obviously  conclude  that  the  process 
of  measuring  predictive  efficiency  of  the  MLR  model  is  meaningless. 
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It  should  be  noted  that  these  results  were  not  unexpected  due  to  the 


low  coefficient  of  determination  (R  = 0.  52)  which  indicated  that  the 
prediction  power  of  the  model  is  low. 


5.  Deviations  of  the  model- -The  deviation  values  of  the  model  are 
shown  in  Table  3-7.  y represents  the  actual  observed  value,  y repre- 


sents the  forecasted  value  and  |y-y  | represents  the  absolute  deviate. 


Table  3-7 

MLR  Model  Deviations 


No. 

y 

y 

1 

5y-y| 

ENGINE 

1 

4 

7.59 

3.59 

2 

3 

8.4 

5.4 

3 

2 

7.86 

5.86 

4 

4 

8.25 

4.25 

5 

6 

8.25 

2.25 

6 

10 

8.43 

1.57 

7 

12 

7,80 

4.20 

8 

11 

7.20 

4.80 

9 

10 

8.25 

1.75 

10 

9 

7.68 

1.32 

11 

14 

7.47 

6.53 

12 

8 

7.35 

0.65 

13 

11 

8. 34 

2.66 

14 

12 

7.53 

4.47 

15 

3 

7.59 

4.59 

WHEEL  ((  BRAKE 

16 

22 

19.71 

2.29 

17 

25 

17. 94 

7.06 

18 

31 

19.53 

11.47 

19 

19 

17.61 

1.39 

20 

23 

18. 18 

4.82 

M 


Table  3-7  (continued) 


No. 

y 

9 

ly-yt 

WHEEL  tc  BRAKE  (continued) 

21 

9 

18.  50 

9.50 

22 

7 

18.67 

11.67 

23 

1 

4.04 

3.04 

COMMUNICATION 

24 

4 

3.97 

0.03 

25 

6 

4.  06 

1.94 

26 

8 

4.  02 

3.98 

27 

9 

4.  11 

4.89 

28 

9 

4.06 

0.06 

29 

13 

12.68 

4.68 

30 

7 

4.08 

3.92 

TOTAL 

124.63 

MAD 

4.  154 

Table  3-8  presents  the  MAD  values  calculated  for  each  class 
using  the  MLR  model 


Table  3-8 

Summary  of  MAD  V alues  - MLR 


CUss 

MAD 

Engine 

3.  592 

Wheel  k Brake 

6.405 

Comnumication 

2. 795 
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The  values  of  Table  3-8  were  subsequently  used  to  compare 
MLR  and  time  series  smoothing  models. 

Time  Series  Results  Presentation  and  Interpretation 

1.  Moving  average- -Table  3-9  presents  the  MAD  results 
that  were  obtained  from  the  computer  output. 

Table  3-9 

Summary  of  MAD  Values  - Moving  Average 


Class  ^ 

N=1 

N=2 

N=4 

N=6 

Engine 

3. 70455 

3. 33750 

3.04000 

Wheel  it  Brake 

3.86956 

4. 11364 

4.32500 

5. 16667 

Communication 

Z. 65217 

2. 72727 

3.22500 

3.61111 

Clas 

N=9 

N=12 

z 

It 

z 

II 

00 

Engine 

2.54815 

2.000 

2. 36296 

Wheel  Si  Brake 

6.46667 

7.  mil 

8.97778 

Communic  a tion 

3.20000 

3.06250 

3.00747 

3.85185 

Conclusions-- 

a.  for  the  engine  class  the  optimal  model  was  found  for 
N^IZ,  when  N is  the  number  of  observations  used  to  compute  the 
moving  average. 

b.  for  the  wheel  and  brake  class  the  optimal  model  was 
found  for  N^l. 
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c.  for  the  communication  components  the  optinnal  model 
was  found  for  N=l. 

An  interesting  conclusion  can  be  drawn  from  the  results 

obtained  for  the  two  last  classes.  When  N=1  the  equivalent  a is 

1 - 1 which  means: 

''1+1' 

a.  No  weight  is  given  to  the  values  of  previous  data  except 
for  the  last  value. 

b.  The  forecasted  demand  for  period  i+1  is  the  dennand  for 

period  i. 

c.  The  interpretation  of  (a)  and  (b)  is  that  the  forecasted 
demand  is  a memoryless  process. 

2.  Exponential  Smoothing --The  minimum  MAD  was  obtained  for  each 
of  the  classes  by  the  remaining  smoothing  models.  This  information 
is  presented  in  Table  3-10. 


Table  3-10 

Minimum  MAD  Values  - Exponential  Smoothing 


Clas 

Type  of 
Smoothing 

a 

Cycle 

Minimum 

MAD 

Engine 

1 

4 

2.51160 

Wheel  k Brake 

1 

1 

2.98078 

Communication 

1 

B 

11 

0. 76464 
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f Conclua ions --The  fact  that  first  order  exponential  smoothing  was 

determined  to  be  optimal  for  all  three  classes  suggests  that  the 
demand  has  no  trend  component  and  the  forecasted  demand  fluctuates 
randomly  relative  to  a certain  constant  that  depends  on  the  specific 
model  (specific  a , and  cycle)  for  each  class. 


Time  Series  Comparison  -Analysis -- 

a.  A comparison  of  minimum  MAO  for  the  optimal  models 
is  given  in  Table  3-11. 


Table  3-11 

MAD  Values  for  Time  Series  Comparison 


^.^''''^’^Optinaal  Model 

Moving  Average 

Exponential 

Smoothing 

Engine 

2.0* 

2.51160 

Wheel  It  Brake 

3.86956 

2. 98078* 

Communication 

2.65217 

0.76474* 

The  results  shown  in  Table  3-11  indicate  that  the  exponential 
smoothing  models  are  "better"  (depends  on  the  predictive  efficiency 
test)  than  the  moving  average  for  the  last  two  classes.  As  for  the 
first  class  (engine)  we  can  conclude  that  the  moving  average  model  is 
more  appropriate;  however,  statistically  significant  conclusions 
cannot  be  drawn  because  of  inappropriate  sample  siaes. 


I 


52 


b.  The  predictive  efficiency  test  of  the  exponential  smooth 
ing  model  for  C.  L.  s 90%  is  shown  in  Table  3-12.  y represents  the 
actual  observed  value,  y represents  the  forecasted  value,  and  y-y 
represents  the  absolute  deviate. 


Table  3-12 

Predictive  Efficiency  Test  of  Exponential 
Smoothing  Model 


Engine 


Wheel  & Brake 


Communication 


Table  3-12  (continued) 


Engine 

Wheel  & Brake 

Communication 

□ 

Time 

n 

9 

A. 

iy-y| 

■ 

a 

9 

|y-9l 

y 

9 

ly-yl 

20 

15 

12.58 

2.42 

♦ 

29 

22.41 

6.59 

* 

B 

8.88 

0.  12 

* 

21 

11 

12.96 

1.96 

* 

24 

26.37 

2.37 

* 

1 

13.23 

1.77 

* 

22 

8 

11.36 

3.36 

* 

30 

24.95 

5.05 

♦ 

B 

6.  14 

0.  14 

* 

23 

12 

12.37 

0.37 

29 

27.98 

1.  02 

* 

8 

8.07 

0.07 

* 

24 

10 

10.32 

0.32 

* 

29 

28.59 

0.41 

* 

9 

8.  15 

0.85 

♦ 

Total  50.  16  65.55  19.1 


MAD  2. 64 

MAD  3.45 

MAD  1.005 

C.I.  90%  = + 5.44 

C.I.  90%  = + 7.  11 

C.I.  90%  = + 2.079 

'*‘Meana  actual  value  "falls"  in  C.I. 

C.L  90%  = K ^„MAD  K ^ = 2.062 

90%  90% 

Conclusions --Table  3-13  summarizes  the  efficiency  levels  obtained 
for  each  class. 


Table  3-13 

Efficiency  Levels  - Summary 


Class 

No.  of  ♦ 

A 

No.  of 

Observations 

B 

Percentage 
A/B  X 100 

m 

Engine 

16 

19 

84.  2% 

90% 

Wheel  k Brake 

17 

19 

89.4% 

90% 

Communic  ation 

16 

19 

84.2% 

90% 

'I 


► 
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From  this  table  we  see  that  the  model  is  most  efficient  for  the  wheel 
and  brake  class.  However,  the  differences  are  too  small  to  conclude 
any  significant  variation  between  the  classes. 

COMPARISON  OF  MLR  MODEL 
AND  TIME  SERIES  MODEL 

1.  Comparison  of  MAD  Values 

Table  3-14  shows  the  different  MAD  values  obtained  for  each 
class  for  the  compared  models. 

Table  3-14 

MAD  Values  - MLR  and  Time  Series  Comparison 


Class  

Model 

MLR 

Time  Series 

Remarks  for 
Time  Series 

Engine 

3.592 

2.0* 

Moving 

Average 

Wheel  It  Brake 

6.405 

2.98078* 

Exponential 

Smoothing 

C omxnunic  ation 

2.785 

0.76474* 

Exponential 

Smoothing 

Conclusion- -The  table  values  indicate  that  in  each  case  the  time  series 
model  is  better  than  the  MLR  model  (using  the  minimum  MAD  as  the 
criterion). 

2.  Comparison  of  efficiency 

The  results  of  the  efficiency  test  of  each  model  disclosed  the 


following! 
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a.  The  efficiency  measurement  for  the  MLR  model  is  mean> 
ingless  (see  page  47,  paragraph  4). 

b.  The  time  series  model  demonstrated  high  predictive 
efficiencies  (84.2%,  89.4%,  84.2%). 

CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  the  results  of  this  chapter  we  conclude  that  the 
time  series  forecasting  technique  is  more  accurate  than  MLR  tech- 
nique. 

The  fact  that  applying  this  technique  is  also  less  complicated 
and  cheaper  than  the  MLR  technique  (less  data,  fewer  variables) 
brings  us  to  the  conclusion  that  time  series  techniques  are  preferable 
to  the  MLR  technique  in  forecasting  future  demand  of  recoverable 
items.  Therefore,  this  forecasting  method  is  used  in  the  next  chapter 
as  an  input  for  development  of  our  general  model. 

It  should  be  noted  again  that  the  conclusions  and  the  recom- 
meiMlation  in  this  chapter  were  based  upon  a small  sample.  The 
methods  and  comparison  techniques  shown  in  this  chapter  can  be  used 
as  an  appropriate  tool  for  a deeper  research  in  this  area. 

I 


1 


CHAPTER  IV 

DEVELOPMENT  OF  THE  MODEL 

Introduction 

As  stated  in  previous  chapters,  the  model  developed  in  this 
research  is  intended  to  determine  two  basic  outputs,  the  quantity  to 
be  recovered  in  the  next  period,  and  the  priority  to  be  assigned  to 
each  item.  These  will  serve  as  inputs  to  the  depot  shop  manager's 
decision  making  process. 


Determination  of  the  Quantity  to  Be  Recovered 

In  determining  the  quantity  of  each  item  to  be  recovered  the 
following  variables  are  considered. 

a.  The  main  input  factor- -the  predicted  consumption  of  each 
item  in  the  next  period  (Cp)  developed  in  Chapter  111. 

b.  The  total  quantity  of  serviceable  items  at  a given  time 

(Q.). 

c.  The  number  of  consumption  periods  (A)  established  by 
Air  Force  policy  for  building  a desired  level  of  serviceable  quantity 
needed  for  future  consumption.  The  future  consultation  is  a function 
of  the  predicted  number  of  war-days. 
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d.  The  number  of  periods  in  which  the  desired  level  of 
serviceable  quantity  will  be  reached  through  the  recovery  process 
(B). 

The  quantity  to  be  recovered  in  the  next  period  is  divided 
into  two  categories: 

a.  The  serviceable  quantity  (Q^)  is  less  than  the  desired 
level  stated  in  the  Air  Force  Policy  (A  * Cp). 

b.  The  serviceable  quantity  (Q^)  is  greater  than  the  desired 

level. 

In  the  first  case,  the  quantity  to  be  recovered  in  the  next  period  (X) 
is  divided  into  two  elements: 

1.  The  quantity  which  compensates  for  the  consumption  in 
the  next  period  (Cp). 

2.  The  relative  part  of  the  desired  serviceable  quantity  level 
which  will  be  built  in  the  next  period 


AC,  - Q, 

B 

in  which: 

A * Cp  is  the  desired  serviceable  quantity  level 
is  the  serviceable  quantity  in  stock 

e 

e e 

(A  * Cp  > Q,)  represents  the  quantity  to  be  repaired  in 
order  to  achieve  the  desired  serviceable 
quantity 
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B is  the  number  of  "building"  periods. 
The  totnl  quantity  to  be  recovered  will  be 


A*  C_ 

+ E. 


Qs 


B 


In  the  second  case  the  serviceable  quantity  is  already  greater  tb*" 
the  desired  level.  In  order  to  prevent  a decreasing  trend  in  the 
serviceable  level,  equation  (2)  was  developed.^ 

X = 2Cp  (2) 

At  first  glance,  it  would  appear  that  in  order  to  maintain  the  existing 
serviceable  quantity,  the  planned  quantity  to  be  recovered  should  be 
equal  to  the  consumption  of  the  next  month  (X  = Cp).  However,  when 
we  consider  the  fact  that  the  actual  recovered  quantity  equals  at  most 
the  planned  quantity  we  can  definitely  conclude  that  using  X s Cp  would 
'**ttlt  in  a continual  deterioration  of  the  serviceable  quantity.  There- 
fore, equation  (2)  is  Used  in  this  case.  The  coefficient  for  Cp  in  this 
aquation  was  found  through  simulation  to  be  the  minimum  satisfactory 
value. 

Determination  of  the  priority  of  each  recoverable  item 
acconqplished  in  accordance  with  I.A.F.  policy  which  states  that  each 

^When  2 ACp  then  X S Cp.  In  addition,  the  quantity  which 
ie  actually  recovered  would  never  be  greater  than  the  planned  quantity 
so  that  the  actual  recovered  quantity  would  not  be  sufficient  to  compen- 
sate for  the  amount  consumed  during  the  next  period.  This  situation 
would  definitely  cause  deterioration  of  the  serviceable  quantity. 
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major  system  is  as  essential  as  any  other.  (The  proposed  model  may 
be  easily  modified  to  reflect  any  changes  in  this  policy. ) 


Two  variables  established  the  priority  for  each  item. 

a.  The  serviceable  quantity  (Q^) 

b.  The  predicted  consumption  in  the  next  period  (Cp). 

The  priority  designator  (Pj)  for  item  j is  defined  as  the  ratio 
between  these  two  variables. 


From  this  equation,  we  see  that  a high  priority  designator  is  a result 
of  a high  Qg  and/or  low  Cp,  and  a low  priority  designator  is  a result 
of  a low  Q,  and/or  high  C^.  A low  priority  designator  therefore 
reflects  a high  essentiality  and  vice  versa. 


(• 


i- 


Figure  4>1 

Flow  Chart  of  tho  Model 
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Note* 


a.  In  considering  the  dispersion  ratio  in  the  model,  we 
made  the  assumption  that  the  minimum  ratio  needed  equals  one.  This 
insures  that  each  user  will  have  at  least  one  item  for  consumption  in 
the  next  period.  This  value  is  totally  subjective  and  is  not  claimed  to 
be  optimal.  Our  intention  in  setting  Dj  s 1 is  to  try  to  prevent  adjust- 
ing between  users  as  a result  of  NORS. 

b.  When  the  calculated  value  of  the  dispersion  ratio  is  less 
than  1,  the  quantity  to  be  recovered  is  changed  to  the  difference 
between  the  number  of  users  and  the  serviceable  quantity.  This  is 
reflected  in  the  following  equations: 


D _ X Q,  _ j 


N 


* e 


X +Q,  = N 

e 

e e 

X » N - 

c.  In  order  to  insure  that  itenos  that  have  a dispersion  ratio 
lower  1 will  be  recovered  in  the  following  period,  a priority 
designator  of  1 is  assigned  to  these  items. 

d.  The  capacity  of  the  depot  shop  was  not  considered  in  the 

model  because  the  model  considers  each  item  separately. 
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The  capacity  of  each  depot  shop  is  taken  into  consideration  while 
comparing  the  total  number  of  hours  needed  to  complete  the  recovery 
program  in  the  next  period.  When  the  total  quantity  of  hours  needed 
is  higher  than  the  depot  shop  capacity,  the  depot  shop  manager  will 
schedule  the  items  according  to  their  priority  designators  (i.  e.  the 
order  of  recovering  will  be  from  the  lower  numbers  to  the  higher 
numbers  until  his  capacity  is  saturated.  However,  when  his  capacity 
is  greater  than  the  total  quantity  of  hours  neede  he  must  finish  all 
of  the  program. 

e.  The  two  factors,  A and  B,  which  reflect  the  desired  level 
of  serviceable  quantity  and  its  "building"  rate,  are  totally  subjective 
and  are  determined  as  stated  in  the  inventory  nnanagement  policy  of 
the  I.  A.  F. 

The  present  policy  states  that  both  A and  B are  equal  to  6,  i.  e.  the 
desired  level  of  serviceable  quantity  is  for  6 months  consumption, 
and  this  level  should  be  achieved  at  the  most  in  six  months  of  the 
recovery  process. 
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CHAPTER  V 

VALIDATION  OF  THE  MODEL 

The  validation  of  the  model  is  conducted  in  two  euccessive 

step*. 

1.  Internal  teste --includes  checking  the  logic  and  consist- 
ency of  the  model. 

2.  External  tests --includes  comparing  model  outputs  with 
real  world  data  through  a computer  simulation.  This  requires  that 
actual  data  be  provided  as  inputs  to  the  model  and  the  resultant  output 
then  be  tested  against  the  criterion  of  being  able  to  achieve  the  speci- 
fied level  of  serviceable  quantity. 

Internal  Test 

The  following  internal  tests  refer  to  the  flow  chart  of  the 
model  in  Chapter  IV  (see  Figure  4-1,  page  6l). 

- Q- 

1.  Step  1 • S ^ 

In  this  aquation  wa  sea  that: 

a.  As  the  predicted  consun^tion  goes  up  (C^)  the  quantity  to 
ba  racovarad  (Xj)  goas  up  also  and  vies  varsa. 
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b.  As  the  serviceable  quantity  (Q^)  goes  up  the  quantity  to 
be  recovered  (Xj)  goes  down,  and  vice  versa. 

c.  As  the  number  of  consumption  periods  (A)  goes  up  the 
quantity  to  be  recovered  also  goes  up  and  vice  versa. 

d.  As  the  "building"  period  (B)  goes  up  the  quantity  to  be 
recovered  (Xj)  goes  down,  and  vice  versa. 


2.  Step  2 


In  this  equation  we  see  that  as  the  serviceable  quantity  (Q^)  goes  up 
the  priority  designator  (Pj)  value  goes  up  also  (i.  e.  low  essentiality 
of  the  item  and  vice  versa.  As  the  predicted  consumption  (C^)  goes 
up  the  priority  designator  (Pj)  value  goes  down  (i.  e.  high  essentiality 
of  item)  and  vice  versa. 

3.  Steps  (3)  (XjSQf)  and  (4)  (X  = Q^)  insure  that  the  quantity 
to  be  recovered  will  never  exceed  the  total  repairable  quantity. 

4.  Steps  (5)  (Dj  > 1),  (6)  (X  s N - Q^)  and  (7)  (Pj  = 1)  insure 
that  whenever  Oj  is  less  than  1,  the  quantity  to  be  recovered  (Xj)  will 
be  such  that  each  user  will  have  a quantity  of  one  of  that  item,  and  it 
will  be  assigned  the  highest  essentiality  level  (Pj  « 1). 

External  Test  j 

The  following  flow  chart  describes  the  external  process  of 
validating  the  model  through  computer  simulation  techniques. 


i 
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Figure  5-1 

Simulation  Flow  Chart 


Explanation* 

1.  The  proceas  starts  with  two  actual  values  of  a certain 
item  Qg- -the  serviceable  quantity,  and  Qj.--the  repairable  quantity 
(see  Steps  (1)  and  (2)). 

2.  In  Step  (3)  the  quantity  to  be  recovered  (Xj)  and  its 
priority  (Pj)  are  determined  by  the  model. 

3.  The  actual  recovered  quantity  (Step  (5))  is  derived  from 
an  experimenl'^l  distribution  that  was  built  from  historical  data 
(description  of  the  distribution  is  given  on  the  following  pages),  and 
influences  the  new  serviceable  quantity  and  repairable  quantity. 

4.  The  actual  consumed  quantity  (Step  (7))  is  derived  from  a 
Poisson  distribution  with  a mean  (X)  that  was  computed  from  histori- 
cal data  (for  more  details  see  page  71),  and  influence  the  predicted 
cons uzx^tion,  the  repairable  quantity  and  the  serviceable  quantity. 


Actual  Recovered  Quantity  Distribution 

This  distribution  gives  the  actual  quantity  repaired  as  a func- 
tion of  the  quantity  which  was  planned  (X^).  It  was  derived  from  500 
double  data  points  (planned  vs  actual)  which  enabled  us  to  build  the 
experimental  distribution. 

The  500  data  points  were  gathered  by  taking  20  months  of 
past  daU  (planned  vs  actual  quantity)  for  25  different  items  from  each 
class.  The  resulto  are  shown  in  Table  5-1  and  the  histograms  pre- 
sented in  Figures  5-2  thru  5-4. 


4 


\ 


o 
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Table  5-1 

Grouped  Data  of  the  i^ctual  Repaired  Quantities 


Quantity 

Class 

** 

X 

X-1 

X-2 

X-3 

X-4 

X-5 

X-6 

X-7 

X-8  and 
more* 

Engine 

252 

88 

59 

37 

28 

14 

8 

9 

5 

Wheel  it 
Brake 

255 

86 

56 

42 

26 

13 

9 

10 

4 

Commun- 

ication 

1 

249 

85 

62 

41 

25 

15 

10 

9 

4 

E 

756 

259 

177 

120 

79 

42 

27 

28 

13 

The  extreme  points  of  these  observations  are  (X-8)  and  (X). 

% 

The  number  of  observations  beyond  the  value  of  (X-8)  is  small  and 
therefore  is  added  to  this  group. 

’*‘*The  depot  never  recovers  more  than  the  planned  quantity. 


Quantity  A:tually  Repaired 


Figure  5>2 

Oiatribution  of  Quantity  Actually  Repaired:  ' 
Engine  Class 
68 


Number  of  Obaervations  Number  of  Observationa 


Figure  5-3 

Diatribution  of  Quantity  Actually  Repaired: 
Wheel  and  Brake  Claas 


From  Table  5-1  and  the  distribution  figures,  we  can  con- 
clude that  there  is  no  significant  difference  in  the  distributions  of 
the  three  classes,  hence  they  can  be  combined  into  a common  dis- 
tribution w^Ach  will  reflect  the  distribution  of  the  actually  recovered 
quantity  as  a function  of  the  planned  quantity  for  the  whole  popula- 
tion, regardless  of  the  different  classes.  The  common  distribution 
is  depicted  in  Figure  5-5. 


Figure  5-5 

Combined  Distribution  of  Quantity  Actually  Repaired 


From  this  histogram  we  can  establish  the  following  distribution 


table.  Table  5-2. 


Table  5-2 

Cumulative  Distribution  Table 


Var. 

X 

X-1 

X-2 

X-3 

X-5 

HI 

X-8 

f(x) 

0.  50 

0.  17 

0. 12 

0.08 

0.05 

0.03 

0.02 

0.02 

0.01 

F(x) 

0.50 

0. 67 

0.79 

0.87 

0.92 

0.95 

0.97 

0.00 

R.N. 

H-SO 

51+67 

68+79 

80+87 

88+92 

93+95 

96+97 

98+99 

100 

The  Actual  Consumption- -Poisson  Distribution 

As  previously  stated,  the  distribution  of  the  actual  consump- 
tion can  be  described  as  a Poisson  process. 

A random  variable  X is  said  to  have  a Poisson  distribution  if 
its  probability  distribution  can  be  written  as 

ki-X 

P|x  = k(  =P^(k)  =-i-—  . 

* k! 

where  X is  a positive  constant  and  k is  any  nonnegative  integer.  It  is 
evident  that  Pjj(k)  is  nonnegative ; and  it  is  easily  shown  that 

jki-k 


The  Poisson  distribution  is  often  used  in  Operation  Research. 
Heuristically  speaking,  this  distribution  is  appropriate  in  many  situ- 
ations where  an  "event"  occurs  over  a period  of  time,  like  the  arrival 
of  a customer.  The  number  of  Customer  Arrivals  in  a fined  time  is 
often  assumed  to  have  a Poisson  distribution.  Sinoilarly  the  demand 
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for  a given  product  is  also  often  assumed  to  have  this  distribution 
( :318).  In  our  case  the  number  of  customer  arrivals  is  the  demand 


and  the  fixed  time  is  a one  month  period. 

As  mentioned  previously,  our  population  of  interest  is 
assumed  to  have  achieved  a stable  trend  of  demand  behaviors.  Thus 
we  can  compute  the  mew  of  the  demand  and  assume  that  it  is  the 
mean  (X.)  of  the  Poisson  distribution. 


CHAPTER  VI 


DESCRIPTION  OF  THE  MODEL  AND  THE 
SIMULATION  PROCESS 

This  chapter  describes  how  the  processes  that  were  devel- 
oped in  Chapters  IV  and  V are  combined  into  one  model. 

Figure  6-1  presents  a flow  chart  that  describes  this  com- 
bined model. 

Explanations 

1.  The  initial  data  that  starts  the  process  is: 

a.  Q^--the  serviceable  quantity  at  this  point  in  time. 

b.  Q, — the  repairable  quantity  at  this  point  in  time. 

c.  Cp--the  forecasted  demand  for  the  next  period  (com- 
puted through  time  series  analysis — See  Chapter  HI). 

d.  N--the  number  of  users. 

e.  Alpha--the  smoothing  coefficient  for  this  item  (see 

Chapter  HI). 

f.  QLamda--the  mean  of  the  demand  distribution 

(Poisson). 

g.  NPNTS--the  number  of  simulation  loops. 
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Figure  6-1 

Simulation  Model  Flow  Chart 


2.  = o • Z + (l-a)Cp.  This  equation  describes  a first 
order  smoothing  function,  2uid  is  common  to  the  three  classes  invest- 
igated. The  difference  between  the  classes  is  reflected  by  different 
Alpha  (or)  values  (for  more  details  see  Chapter  III,  page  51). 

3.  The  random  number  (R.N.)  for  the  actually  repaired 
quantity  is  generated  by  using  the  random  number  generator  on 
CREATE  (ULIB)  GRADLIB/TSS,  R.  By  using  Seed  = +1  we  are  able 
to  generate  the  same  sequence  of  R.N.  for  all  the  cases  investigated 
(7:Chap  7). 

4.  The  actually  repaired  quantity  is  chosen  according  to  the 
R.N.  and  the  distribution  which  is  described  in  Chapter  V,  Table  No. 
5-2. 

5 and  6.  The  actually  consumed  quantity  is  generated  using 
the  Poisson  distribution  generator  which  exists  on  CREATE. 
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CHAPTER  Vn 

SIMULATION  RESULTS  AND  CONCLUSIONS 

Introduction 

In  this  chapter  we  describe  the  results  of  the  simulation. 

! 

) 

^ The  process  was  simulated  for  each  class  of  items  (engine, 

wheel  and  brake,  conununication).  In  each  class  we  examined  six 
j different  initial  situations,  which  can  be  classified  into  two  main 

groups: 

$ 

^ a.  high  serviceable  quantity  and  low  repairable  quantity 

^ b.  low  serviceable  quantity  and  high  repairable  quantity. 

i 

The  first  group  is  a test  case  examining  the  ability  of  the  model  to 
prevent  deterioration  of  the  existing  level  of  serviceable  quantity  (see 
Chapter  IV,  page  59). 

The  second  group  is  a test  case  which  examines  the  ability 
of  the  model  to  build  a serviceable  quantity  up  to  a desired  level  (A  • Q) 
within  a desired  number  of  periods  (B).  Each  group  is  divided  into 
three  situations: 

a.  K < Cp,  which  means  that  the  mean  of  the  future  consump- 
tion distribution  is  lower  than  our  initial  prediction.  | 

I 

4 

I 

5 * 

L 
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b.  X.  = Cp,  which  meains  that  the  mean  of  the  future  consump- 
tion distribution  equals  our  initial  prediction  (which  is  the  case  of 
reaching  a stable  trend  of  behavior). 

c.  X > Cp  which  means  that  the  mean  of  the  future  consump- 
tion distribution  is  higher  than  our  initial  prediction. 

The  differentiation  between  the  three  classes  of  items  is 
done  by  using  the  three  different  smoothing  factors  (a's)  which  were 
obtained  in  the  time  series  analysis  (see  Chapter  III,  page  51).  Those 
three  a's  lead  to  the  difference  between  the  predicted  future  consump- 
tion (Cp)  for  each  class. 

Cp  = Of  • Z + (l-o()Cp 


Initial  Simulation  Data 

Table  7-1  gives  the  initial  data  points  used. 
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Table  7-1 

Initial  Data  Points  of  the  Simulation 


Class 


Engine 


Wheel  & 
Brake 


Commun-  | 50 
ication 


NPNTS  IS  the  number  of  simulation  periods,  which  in  our  case  is  100 
months  (more  than  8 years  into  the  future).  The  high  number  was 
chosen  in  order  to  prevent  the  influence  of  initial  conditions  and  to 
study  the  behavior  of  the  model  in  the  long  run. 

A = B = 6 is  chosen  because  the  existing  policy  of  the  Air  Force  is 
to  build  a serviceable  quantity  for  six  months  of  consumption  in  a 
period  of  six  months. 


In  addition  to  the  previous  situations  which  were  examined  in 
accordance  with  the  existing  Air  Force  policy,  we  checked  another 


situation  with  different  stated  policies: 
a.  A • 6,  B * 3. 
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This  policy  requires  building  a serviceable  quantity  for  six  months 
of  consumption  in  a period  of  3 months  (a  higher  rate  of  building  the 
desired  level  in  comparison  to  the  existing  Air  Force  policy), 
b.  A = 12,  B = 6. 

This  policy  requires  building  a serviceable  quantity  for  12  months  of 
consumption  in  a period  of  6 months  (a  higher  desired  serviceable 
level  in  comparison  to  the  existing  Air  Force  policy). 

The  simulation  computer  program  auid  the  simulation  results 
(quantitative  and  descriptive-graphs)  are  presented  in  Appendix  A. 

Results  Analysis 

The  results  analysis  refers  to  the  initial  data  which  is  pre- 
sented in  Table  7-1.  From  the  results  we  can  see  that  the  differences 
between  the  three  groups  under  investigation  are  minor  and  do  not 
influence  the  model's  ability  to  operate  in  the  expected  manner.  This 
fact  enables  us  to  analyze  the  results  of  one  group  and  make  a general 
conclusion  about  the  behavior  of  the  model  concerning  all  three 
classes. 

Situation  #1 

Q,  = 50 
Q,  = 5 

S'  * 

\ s 2 
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or  = 0.  5 


NPNTS  = 100 
A = 6 
B = 6 
N = 5 

This  initial  situation  represents  a high  level  of  serviceable 
quantity,  a low  level  of  repairable  quantity  and  a low  average  con- 
sumption level  relative  to  the  initial  Cp.  As  expected,  we  see  that 
the  level  of  serviceable  quantity  is  not  significantly  changed  in  com- 
parison to  the  initial  condition.  This  allows  us  to  assume  that  the 
model  does  not  "permit”  deterioration  of  the  high  serviceable  quantity. 

Situation  #2 

Qs=2 

Qy  = 30 

Cp  = 5 

X =2 
a =0.5 
NPNTS  = 100 
A = 6 
B = 6 
N = 5 

i 

This  initial  situation  represents  a low  level  of  serviceable  | 

quantity,  a high  level  of  repairable  quantity  and  a low  average  j'  ■ 
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consumption  level  relative  to  the  initial  Cp.  As  expected  we  see  that 
the  model  "builds"  the  serviceable  quantity  to  the  desired  level 
(Qg  = 6Cp)  in  six  months,  and  then  continues  to  improve  this  quantity 
as  much  as  possible  (the  upper  limit  is  the  total  quantity  of  this  item). 


Situation  #3 


Qg  = 50 
Q,  = 5 


Cp=  5 
X = 5 


a =0.5 


NPNTS  = 100 


A = 6 


B = 6 


N = 5 


This  initial  situation  represents  a high  serviceable  quantity, 
a low  repairable  quantity  and  a median  level  of  consumption  equal  to 
the  initial  Cp.  As  expected,  we  see  that  the  model  behaves  very  simi> 
larly  to  Situation  No.  1. 


Q,  » 2 


Q,  ■ 30 


• 5 
P 


a 


r 


r I 


X = 5 


a s 0.  5 


NPNTS  « 100 


A s 6 


B > 6 


N « 5 


Thia  initial  situation  represents  a low  serviceable  quantity, 
a high  repairable  quantity  and  a median  level  of  consumption  equal  to 
the  initial  Cp.  As  expected,  we  see  that  the  model  behaves  very  simi- 
larly to  Situation  #2. 


Situation  #5 

Q,  « 50 
Qr«  5 
Cp*  5 
X * 8 
a B 0. 5 
NPNTS  ■ 100 
A ■ 6 
B ■ 6 
N B 5 

This  initial  situation  represents  a high  serviceable  quantity, 
a low  repairable  quantity  and  a high  average  consumption  level  relative 


to  the  initial  Cp.  A*  expected,  we  eee  that  the  model  behaves  very 
einularly  to  Situations  #1  and  3. 


r 


J 


Situation  #6 

Q.  = 2 
Q,  = 30 
Cp=  5 
X = 8 
a =0.5 
NPNTS  = 100 
A = e 

B = 6 
N = 5 

This  situation  represents  a low  serviceable  quantity,  a high 
repairable  quantity  and  a high  average  consumption  level  relative  to 
the  initial  Cp.  As  expected,  we  see  that  the  model  behaves  very  simi- 
larly to  Situations  #2  and  4. 

Situation  #7 

Qf  « 2 

Q,  - 30 

S'  * 

X ■ 2 
« ■ 0.5 


iM 


I 

} 

I 

5 


NPNTS  « 100 
A = 12 
B s 6 
N = 5 

This  situation  is  similar  to  Situation  #2  in  its  initial  condi- 
tions accept  for  the  change  in  the  value  of  A.  We  examined  this  situ- 
ation to  see  how  a change  in  the  desired  level  of  serviceable  quantity 
(A  * Cp)  is  reflected  in  the  behavior  of  the  model.  According  to  the 
results  we  see  that  the  rate  of  building  the  desired  serviceable  quant- 
ity is  higher  than  the  case  in  Situation  #2;  i.  e. , n this  situation  we  see 
that  after  a six  month  period  we  reach  a serviceable  quantity  of 
Qg  * 27.4  while  in  Situation  #2  we  reached  a quantity  of  Qg  = 14.  5. 
However,  it  should  be  stated  that  in  the  long  run  the  process  stabilizes 
for  the  two  situations  around  the  same  level  of  serviceable  quantity 

(Qgat  28). 


r 


B = 3 


N = 5 

This  situation  is  similar  to  Situation  #2  in  the  initial  condi- 
tions, except  for  the  change  in  the  value  of  B.  We  examined  this 
situation  to  see  how  the  model  responds  to  a change  in  the  desired 
"building  period"  (B).  According  to  the  results,  we  see  that  the  rate 
of  building  the  serviceable  quantity  increases;  i.  e. , in  this  situation 
it  was  observed  that  after  a three  month  period  we  reach  a service- 
able quantity  of  Qg  s 18.  3 while  in  Situation  #2  we  reached  a quantity 
of  Qg  s 11.7.  However,  it  should  be  stated  that  in  the  long  run  the 
process  stabilises  for  both  situations  around  the  same  level  of  service- 
able quantity  (Qg  a 28). 

General  Conclusions 

According  to  the  results  analysis  we  see  that  the  model 
behaves  as  expected. 

1.  In  all  cases  where  we  specified  initial  conditions  of  high 
serviceable  quantities  and  low  repairable  quantities  the  end  result 
was  that  the  serviceable  qusAtity  was  not  deteriorating  and,  moreover, 
the  maximum  possible  serviceable  quantity  was  achieved. 

2.  In  all  cases  where  we  had  initial  conditions  of  low  service- 
able quantities  and  high  repairable  quantities,  the  model  built  the  , 
serviceable  quantity  and  reached  the  desired  level  within  the  specified 

amount  of  time,  (i.  e.  according  to  the  values  of  A and  B).  In  the  long 
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run,  the  model  continued  to  build  the  serviceable  quantity  and  stabi- 
lised it  around  the  maximum  possible  level. 

3.  In  establishing  priority  designators  for  all  the  cases  we 
were  able  to  Identify  a similar  trend  of  behavior.  When  the  service- 
able quantity  (Q,)  goes  up,  the  repairable  quantity  (Q^)  goes  down  and 
the  priority  value  goes  up,  which  means  a lower  essentiality  is 
assigned  to  that  particular  item. 

The  priority  designators  give  the  depot  shop  manager  an 
indication  as  to  how  he  should  schedule,  within  the  constraint  imposed 
by  shop  capacity,  the  repair  of  various  items.  We  can  say  that  in  the 
case  of  building  serviceable  quantity  (low  Qg,  high  Q, ) the  priority 
value  influences  the  scheduling  of  the  item  in  the  shop's  workload. 
According  to  the  results  we  see  that  in  the  short  run  we  have  low 
priority  values  (high  essentiality)  in  successive  periods.  This  means 
that  this  item  will  be  scheduled  for  repair  in  successive  periods, 
while  in  the  long  run  we  see  that  the  number  of  successive  periods  in 
which  low  priority  items  are  repaired  goes  down  (indicating  less  suc- 
cessive scheduling  of  repair  of  the  item  in  the  shop  workload). 

In  the  case  of  initial  high  serviceable  quantity  and  low  repair- 
able quantity  we  see  that  we  obtain  a high  priority  designator  which 
indicates  that  the  item  is  not  expected  to  be  scheduled  in  the  shop 
(dapanda  on  tha  capacity  of  tha  shop). 
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According  to  the  conclusions  stated  above  we  can  say  that  the 
model  which  was  developed  in  this  thesis  responds  as  expected  for  the 
different  situations  that  were  investigated.  Thus,  the  model  is  a 
viable  tool  for  implementation  in  the  depot  recovery  process  for 
repairable  items  in  the  Air  Force. 

Recommendations  for  Further  Research 

The  following  areas  are  recommended  for  further  research: 

a.  An  investigation  of  the  sensitivity  of  the  model  to  differ- 
ent dispersion  ratios. 

b.  An  investigation  of  the  influence  of  the  distribution  of 
actually  repaired  quantities  on  the  "building"  capability  of  the  model 
and  the  sensitivity  of  the  overall  model  to  this  distribution. 

c.  Constructicm  of  a similar  model  incorporating  non-linear 
relationships  between  the  variables. 
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APPENDIX  A 

SIMULATION  PROGRAM  AND  RESULTS 
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non  to  • 

SOON  •= ( UL I n )CRADL 1 n/TSS,R 

000 15C 

THIS  PROCRAM  COMPUTES  THE  RESIIITS 

FOR  THE 

ooo?ac 

simulatiom  process  or  the  model 

WHICH 

0 0023C 

DETERMINES  THE  QUANTITY  TO  BE  RECOVERED 

00030C 

IN  THE  DEPOT  IN  THE  NEXT  PERIOD 

AND  ITS 

00035C 

PRIORITY 

■ 

0 0 0-10 

DIMENSION  PERtOD(100).OSl(100> 

,OR1(100) 

000-«5 

DIMENSION  Xl(100),Yl(100),Pl(ln0) 

00050 

DIMENSION  ZKIOOI.CPKIOO) 

* 

00055 

REAL  NPNTS 

00060 

READ, OS, QR. CP. alpha, QLAHDA. 

NPNTS 

00065 

XOS=OS 

00070 

XORsOR 

00075 

XCPsCP 

OOOOO 

Y = 0 

00065 

z = o 

00090  ' 

Na5 

00095 

PRINT  7001 

0 n * 0 

7001  FORMAT! //////13X, '•*••**  SIMULATION 

«6X,"0S",5X,"QR*',6X,"X",6X," 

P", 

6X, " Y" 

t6X."Z",6X,«CP'',//) 

20 

DO  1 J^l, NPNTS 

00125 

1 PERIOD(J>aJ 

. 00130 

DO  73  1=1, NPNTS 

00135 

OSaOS* Y-Z 

00110 

ORsOR*Z-y 

00145 

iriOS.LT.O. )G0  TO  2 

00190 

IF<OR.LT.O, >00  TO  3 

00155 

60  TO  4 

00160 

7 OS=0. 

00165 

00  TO  4 

00170 

3 &RaO. 

00175 

4 P=OS/CP 

. 00180 

E=6*CP 

00185 

ir(OS.CE.E)CO  TO  5 

00190 

X=CP*(6*CP-0S)/6 

00195 

GO  TO  6 

00700 

5 X=2*CP 

00205 

6 IFtX.LT. 1.0)00  TO  10 

) 

00710 

7 IF(X.CE.QR>CO  TO  8 

00715 

D>( X«OS)/N 

00770 

IFCII.LE.l  ,0)00  TO  9 

00775 

Cn  TO  It 

00730 

8 X*0R 

* • .'r** 

90 

1 

00235 

CO  TO  11 

0 0 2 4 0 

9 

X = N-OS 

00245 

P=1 . 0 

00250 

00  TO  11 

00255 

10 

X=1.0 

00260 

11 

PRINT  7002, OS, OR. 

X,  P 

, Y,7,CP 

00265 

OSKDaOS 

00270 

ORl ( I ) sOR 

00275 

XI ( I lax 

00280 

PI ( I )»P 

00285 

Yl( 1 )sY 

00290 

Zl( 1 >>i 

00295 

CPK  I )aCP 

00300 

7002 

rORHAT(5X.7(r4,l, 

3X)/) 

00305 

SEEDS  ♦l.o 

00310 

RV  > RND( SEED)*100»1 

00315 

IF(RV.GT,50,0)GO 

TO 

12 

00320 

YaX 

00325 

GO  TO  9 

00330 

12 

irCRV. 07.67.0)00 

TO 

13 

00335 

Y=X-1 

00340 

GO  TO  9 

00345 

13 

IF(RV.Or.79.0)GO 

TO 

14 

00350 

YaX-2 

00355 

GO  TO  9 

- 

00360 

14 

IFCRV.OT. 87.0)00 

TO 

15 

00365 

YaX-3 

00370 

GO  TO  9 

00375 

15 

lF(RV.Or.92.0)00 

TO 

16 

00380 

YaX-4 

00385 

GO  TO  9 

00390 

16 

lF(RV.Or.95.0)GO 

TO 

17 

00395 

YaX-5 

00400 

60  TO  9 

00405 

17 

IF(RV.GT.97.0)00 

TO 

18 

0 0410 

Y«X-6 

00415 

60  TO  9 

00420 

18 

IFIRV.OT. 99.0)00 

TO 

19 

00425 

Y-X-7 

00430 

60  TO  9 

00435 

19 

Y«X-8 

00440 

20 

CONTINUE 

00445 

IF( Y.LE.P. )00  TO 

21 

00450 

Y»Y 

00455 

GO  TO  22 

0 0 460 

21 

Y = 0 

0 0 46^ 

22 

Z=  POISSON  (QLAMDA,!.) 

0 0 4 7 0 

CP  = ALPHA»2*(1-A1,PHA)*CP 

00475 

23 

CONTINUE 

00  480 

CALL 

11 S T A R T 

0 0485 

CALL 

USET( "SMALL"  ) 

0 0490 

CALL 

IIDlMFN(12.,12.,"AriT/LSG 

0 0495 

OTOL l 3AN0  BRAUN\"  ) 

00500 

CALL 

UDAREA(0.,0.5,0.,11.) 

00505 

call 

UOUTLN 

00510 

CALL 

UDAREA(1.,7.5,1,,4.25) 

00515 

call 

USET( "XROTH"  ) 

00520 

CALL 

USET( "YBOTH"  ) 

005?5 

CALL 

UPSET("XLAB",.  "PERIODV") 

00530 

CALL 

UPSET< " YLAB"  , "P\"  ) 

00535 

CALL 

UPLOTl (PERIOD, PI, NPNTS) 

00540 

CALL 

LAB(1.5,4.25,l..55,,"P\" 

00545 

CALL 

UDAREA(1.,7,5,5,75,9,0) 

00550 

CALL 

UPSET("YLA8", "OS  OR\") 

00555 

CALL 

UPLOTl(PERIOD,OSl,MPNTS) 

00560 

CALL 

UPSET("SETDASH",34. » 

00565 

CALL 

UHOVE(PERIOD(l»,ORl<l)> 

00570 

DO  24 

J=l, NPNTS 

00575 

24 

CALL 

UPENl(PER10n(J),0Rl(J)," 

00580 

CALL 

LAB(l,5,9,n,l.,55.,"0S\" 

00555 

CALL 

LAB<1,5,9.0,3.,34.,"OR\" 

00590 

‘CALL 

USET< "DEVICE"  ) 

00595 

CALL 

UM0VE(6. ,10 . ) 

00600 

CALL 

UDO!T("LF01") 

00605 

CALL 

UPRNTl("OSa","TEXT" ) 

#0610 

CALL 

UPRHT1(X0S,"REAL") 

00615 

CALL 

UMOVE(6.,10.  ) 

00620 

CALL 

UDO!T{"LF02" ) 

00625 

CALL 

UPRNTK  "ORa",  "TEXT"  ) 

00630 

CALL 

UPRNTK  XOR,  "REAL"  ) 

00635 

CALL 

UHOyE(6..10.  ) 

00640 

CALL 

UOOIT<"Lr03") 

00645 

call 

UPRMT1("CP»", "TEXT" ) 

^0650 

CALL 

UPRHTIC  XCP, "REAL" ) 

00655 

CALL 

UHOVE(6.,10. > 

00660 

CALL 

UOOI T< "Lr04" ) 

00665 

CALL 

UPRNTK  "alphas".  "TEXT") 

60670 

CALL 

UPRNTK  alpha,  "REAL"  ) 

00675 

CALL 

UHOVE(6,.10. ) 

00680 

CALL 

UDOl T( "LFn5" ) 

770 


) 


DASH"  ) 
) 

) 
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0 0 6 M 5 

CALL 

UPRNTl("OLAHnA  = ",**TEXT"> 

00600 

CALL 

UPRNTKOLAMnA.-REAL") 

00693 

CALL 

UEND 

00700 

STOP 

00703 

END 

00710 

SUBROUTINE  LABI  X, Y, A 1 , OSH, LABI ) 

00713 

HsAt 

00720 

CALL 

USET(  "DEV  1 Ce'"  ) 

00725 

CALL 

UHOVECX, Y) 

00730 

1 CALL 

UD01T("UP0l" ) 

00733 

HaN-l 

00740 

iriN. 

GT. 0)00  TO  1 

00745 

CALL 

UPSET(-SETDASH-,DSH) 

00750 

CALL 

USETI-DASH" ) 

00755 

call 

UPEN1(1,.0,, "RELATIVE") 

00760 

CALL 

UDOI T< "SPOZ" ) 

00765 

CALL 

UPRNT1(LAD1,"tEXT") 

00770 

CALL 

USET(-LINE"), 

00775 

CALL 

UPSET("S£TDASH",56, ) 

00700 

CALL 

USET<"VIRTUAL") 

00705 

RETURN 

00790 

END 
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*****  SIIIULATIOII  RESULTS  ***** 


INITIAL  CONDITIONS 

QS-50.0  QR-  S.O  CP-  S.O  ALPHA-0.00  LAtlOA-  2.0 


OS 

QR 

X 

P 

Y 

Z 

CP 

50.0 

5.0 

5.0 

10.0 

0. 

0. 

5.0 

55. 0 

0. 

0. 

55.0 

5.0 

0. 

1.0 

52.0 

3.0 

3.0 

20.0 

0. 

3.0 

2.6 

50.0 

5.0 

5.0 

11.1 

3.0 

5.0 

4.5 

49.0 

6.0 

6.0 

11.9 

3.0 

4.0 

4.1 

53.0 

2.0 

2.0 

21.9 

6.0 

2.0 

2.4 

55.0 

0. 

1.0 

13.6 

2.0 

0. 

0.5 

53.0 

2.0 

2.0 

21.2 

1.0 

3.0 

2.5 

50.0 

5.0 

5.0 

13.5 

1.0 

4.0 

3.7 

55.0 

0. 

0. 

74.3 

5.0 

0. 

0.7 

53.0 

2.0 

2.0 

30.3 

0. 

2.0 

1.7 

52.0 

3.0 

2.3 

45.2 

0. 

1.0 

1.1 

51.0 

4.0 

2.1 

49.5 

0. 

1.0 

1.0 

49.  1 

5.9 

5.9 

14.4 

2.1 

4.0 

3.4 

47.0 

8.0 

7.8 

12.  1 

1.9 

4.0 

3.9 

49.8 

5.2 

5.2 

12.5 

6.8 

4.0 

4.0 

50.0 

5.0 

4.8 

20.9 

2.2 

2.0 

2.4 

52.8 

2.2 

2.2 

41.3 

3.8 

1.0 

1.3 

30.8 

4.2 

3.7 

27.4 

0. 

2.0 

1.9 

51.5 

3.5 

3.5 

18.6 

3.7 

3.0 

2.8 

52.0 

3.0 

1.1 

93.8 

0.5 

0. 

0.6 

51.0 

4.0 

1.8 

56.0 

0. 

1.0 

0.9 

50.8 

4.2 

3.6 

28.5 

1.8 

2.0 

1.8 

52.4 

2.6 

2.3 

45.3 

2.6 

1.0 

1.2 

52.7 

2.3 

2.  1 

51.1 

1.3 

1.0 

1.0 

49.8 

5.2 

5.2 

19.  1 

0.  1 

3.0 

2.6 

54.0 

1.0 

1.0 

40.9 

5.2 

1.0 

1.3 

50.0 

5.0 

5.0 

14.4 

0. 

4.0 

3.5 

51.0 

4.0 

4.0 

16.5 

4.0 

3.0 

3. 1 

55.0 

0. 

0. 

88. 9 

4.0 

0. 

0.6 

54.0 

1.0 

1.0 

58.5 

0. 

1.0 

0.9 

54.0 

1.0 

1.0 

92.3 

0. 

0. 

0.2 

53.0 

2.0 

1.7 

63.3 

0. 

1.0 

C.8 

53.0 

2.0 

1.0 

16.6 

0. 

0. 

0.2 

52.0 

3.0 

1.7 

62.4 

0. 

1.0 

0.8 

46.7 

8.3 

8.3 

9.4 

0.7 

6.0 

5.0 

54.0 

1.0 

1.0 

30.  1 

8.3 

1.0 

1.8 

52.0 

3.0 

3.0 

26.5 

0. 

2.0 

2.0 

51.0 

4.0 

4.0 

14.2 

3.0 

4.0 

3.6 

54.0 

1.0 

1.0 

35.6 

4.0 

1.0 

1.5 
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‘iS 

gii 

V 

4b 

P 

Y 

Z 

CP 

54. 0 

1.0 

1.0 

48.9 

1.0 

1.0 

1.1 

50.0 

5.0 

5.0 

14.6 

0. 

4.0 

3.4 

53.0 

■Sa 

2.0 

35.  7 

4.0 

1.0 

1.5 

51.0 

mjm 

4.0 

14.6 

2.0 

4.0 

3.5 

55.0 

BsH 

0. 

78.6 

4.0 

0. 

0.7 

53.3 

2.0 

30.5 

0. 

2.0 

1.7 

52.0 

3.0 

2.3 

45.3 

0. 

1.0 

1.1 

53.3 

1.7 

1.7 

51.8 

2.3 

1.0 

1.0 

53.3 

1.7 

1.0 

58.8 

0. 

0. 

0.2 

52.3 

2.7 

1.7 

62.2 

0. 

1.0 

0.8 

50.0 

5.0 

5.0 

14.8 

1.7 

4.0 

3.4 

48.0 

7.0 

6.  1 

15.6 

1.0 

3.0 

3.  1 

51.1 

3.9 

3.9 

23.  1 

5.  1 

2.0 

2.2 

51.1 

3.9 

1.0 

15.5 

0. 

0. 

0.4 

52.  1 

2.9 

1.0 

68.6 

1.0 

0. 

0.  1 

52.1 

2.9 

1.6 

63.8 

1.0 

1.0 

0.8 

52.8 

2.2 

1.9 

54.3 

1.6 

1.0 

1.0 

51.7 

3.3 

3.3 

19.9 

1.9 

3.0 

2.6 

54.0 

1.0 

1.0 

41.0 

3.3 

1.0 

1.3 

54.0 

1.0 

1.0 

04.8 

0. 

0. 

0.3 

52.0 

3.0 

3.0 

31.5 

0. 

2.0 

1.7 

49.0 

6.0 

13.9 

1.0 

4.0 

3.5 

50.0 

5.0 

4.6 

21.7 

3.0 

2.0 

2.3 

52.6 

2.4 

2.4 

25.5 

4.6 

2.0 

2.  1 

50.6 

4.4 

4.0 

25.2 

0. 

2.0 

2.0 

51.6 

3.4 

2.4 

42.9 

2.0 

1.0 

1.2 

53.0 

2.0 

2.0 

51.0 

2.4 

1.0 

1.0 

55.0 

0. 

1.0 

64.3 

2.0 

0. 

0.2 

55.0 

0. 

0. 

65.4 

1.0 

1.0 

0.8 

53.0 

2.0 

. 2.0 

30.0 

0. 

2.0 

1.8 

54.0 

1.0 

1.0 

46.8 

2.0 

1.0 

1.2 

51.0 

4.0 

4.0 

19.4 

0. 

3.0 

2.6 

52.0 

3.0 

2.7 

39.2 

2.0 

1.0 

1.3 

51.7 

3.3 

3.3 

19.4 

2.7 

3.0 

2.7 

48.7 

6.3 

5.9 

16.6 

0. 

3.0 

2.9 

52.5 

2.5 

2.5 

24.0 

5.9 

2.0 

2.2 

49.5 

5.5 

5.5 

17.5 

0. 

3.0 

2.8 

52.0 

3.0 

3.0 

17.5 

5.5 

3.0 

3.0 

51.0 

4.0 

2.8 

36.6 

0. 

1.0 

1.4 

48.8 

6.2 

6.2 

11.4 

2.8 

5.0 

4.3 

54.0 

1.0 

1.0 

32.6 

6.2 

1.0 

1.7 

54.0 

1.0 

1.0 

63.  1 

0. 

0. 

0.3 

52.0 

3.0 

3.0 

31.2 

0. 

2.0 

1.7 

53.0 

2.0 

2.0 

46.8 

2.0 

1.0 

1.1 

54.0 

1.0 

1.0 

52.6 

2.0 

1.0 

1.0 

54.0 

1.0 

1.0 

53.7 

1.0 

1.0 

1.0 

54.0 

1.0 

1.0 

53.9 

1.0 

1.0 

1.0 

53.0 

2.0 

2.0 

53.0 

0. 

1.0 

1.0 
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*****  SIMULATIOll  RESULTS  ***** 

INITIAL  COMDITIONS 

QS-  2.0  QR-30.0  CP-  5.0  ALPHA-0. 30  LAUDA-  2.0 


ys 

ya 

X 

P 

Y 

Z 

CP 

2.0 

30.0 

9.7 

0.4 

0. 

0. 

5.0 

11.7 

20.  3 

2.0 

11.7 

9.7 

0. 

1.0 

10.7 

21.3 

3.4 

4.  1 

2.0 

3.0 

2.6 

9.  1 

22.9 

7.5 

2.0 

3.4 

5.0 

4.5 

10.6 

21.4 

6.4 

2.6 

5.5 

4.0 

4.1 

IS.  1 

16.9 

4.8 

6.2 

6.4 

2.0 

2.4 

19.9 

12.  1 

1.0 

41.1 

4.8 

0. 

0.5 

17.9 

14.  1 

5.0 

7.2 

1.0 

3.0 

2.5 

17.9 

14.  1 

4.4 

4.8 

4.0 

4.0 

3.7 

22.3 

9.7 

1.5 

30.  1 

4.4 

0. 

0.7 

21.8 

10.2 

3.5 

12.5 

1.5 

2.0 

1.7 

22.3 

9.7 

2.3 

19.4 

1.5 

1.0 

1.1 

21.3 

10.7 

2.  1 

20.7 

0. 

1.0 

1.0 

19.3 

12.7 

3.6 

5.7 

2.  1 

4.0 

3.4 

13.3 

16.7 

5.2 

4.0 

0. 

4.0 

3.9 

15.5 

16.5 

5.4 

3.9 

4.2 

4.0 

4.0 

15.9 

16.  1 

4.8 

6.6 

2.4 

2.0 

2.4 

18.7 

13.3 

2.0 

14.6 

3.8 

1.0 

1.3 

16.  7 

15.3 

3.7 

9.0 

0. 

2.0 

1.9 

17.4 

14.6 

5.5 

6.3 

' 3.7 

3.0 

2.3 

20.0 

12.0 

1.1 

36.0 

2.5 

0. 

0.6 

19.0 

13.0 

1.8 

20.8 

0. 

1.0 

0.9 

18.8 

13.2 

3.6 

10.5 

1.8 

2.0 

1.8 

20.3 

11.7 

2.3 

17.6 

2.6 

1.0 

1.2 

20.7 

11.3 

2.1 

20.0 

1.3 

1.0 

1.0 

17.7 

14.3 

5.2 

6.8 

0.  1 

3.0 

2.6 

21.9 

10.  1 

2.6 

16.6 

5.2 

1.0 

1.3 

17.9 

14.1 

3.9 

5.2 

0. 

4.0 

3.5 

17.9 

14.1 

3.2 

5.8 

2.9 

3.0 

3.1 

21.1 

10.9 

1.2 

34.1 

3.2 

0. 

0.6 

20.1 

11.9 

1.8 

21.7 

0. 

1.0 

0.9 

20.9 

11.1 

1.0 

13.3 

0.8 

0. 

0.2 

19.9 

12.1 

1.7 

23.6 

0. 

1.0 

0.8 

19.9 

12.  1 

1.0 

19.0 

0. 

0. 

0.2 

18.9 

13. 1 

1.7 

22.7 

0. 

1.0 

0.8 

13.6 

18.4 

7.7 

2.7 

0.7 

6.0 

5.0 

20.3 

11.7 

3.6 

11.3 

7.7 

1.0 

1.8 

18.3 

13.7 

3.9 

9.3 

0. 

2.0 

2.0 

18.2 

13.8 

4.2 

5.1 

3.9 

. 4.0 

3.6 

21.3 

10.7 

3.0 

14.0 

4.2  • 

1.0 

1.5 

QS 

QR 

A 

P 

Y 

Z 

CP 

23.4 

3.6 

2.2 

21.2 

3.0 

1.0 

1.1 

19.4 

12.6 

3.6 

5.7 

0. 

4.0 

3.4 

21.0 

11.0 

3.0 

14.1 

2.6 

1.0 

1.5 

19.9 

12.  1 

3.7 

5.7 

3.0 

4.0 

3.5 

23.6 

8.4 

1.4 

33.8 

3.7 

0. 

0.7 

23.0 

9.0 

3.5 

13.2 

1.4 

2.0 

1.7 

23. S 

8.5 

2.3 

20.5 

1.5 

1.0 

1.1 

24.8 

7.2 

2.1 

24.1 

2.3 

1.0 

1.0 

24.8 

7.2 

1.0 

20.4 

0. 

0. 

0.2 

23.8 

8.2 

1.7 

28.3 

0. 

1.0 

0.8 

21. S 

10.5 

6.7 

6.4 

1.7 

4.0 

3.4 

21.2 

10.8 

6.1 

6.9 

2.7 

3.0 

3.  1 

24.4 

7.6 

4.4 

11.0 

5.1 

2.0 

2.2 

24.8 

7.2 

1.0 

56.0 

0.4 

0. 

0.4 

25. 8 

6.2 

1.0 

91.1 

1.0 

0. 

0.  1 

25.8 

6.2 

1.6 

31.5 

1.0 

1.0 

0.8 

26.4 

5.6 

1.9 

27.4 

1.6 

1.0 

1.0 

25.3 

6.7 

5.2 

9.8 

1.9 

3.0 

2.6 

29.5 

2.5 

2.5 

22.4 

5.2 

1.0 

1.3 

30.0 

2.0 

1.0 

13.8 

0.5 

0. 

0.3 

28.0 

4.0 

3.3 

16.9 

0. 

2.0 

1.7 

25.3 

6.7 

6.7 

7.2 

1.3 

4.0 

3.5 

27.0 

5.0 

4.6 

11.7 

3.7 

2.0 

2.3 

29.6 

2.4 

2.4 

14.4 

4.6 

2.0 

2.1 

27.6 

4.4 

4.0 

13.7 

0. 

2.0 

2.0 

28.6 

3.4 

2.4 

23.8 

2.0 

1.0 

1.2 

30.0 

2.0 

2.0 

28.9 

2.4 

1.0 

1.0 

32.0 

0. 

1.0 

53.8 

2.0 

0. 

0.2 

32.0 

0. 

0. 

38.0 

1.0 

1.0 

0.8 

30.0 

2.0 

2.0 

17.0 

0. 

2.0 

1.8 

31.0 

1.0 

1.0 

26.9 

2.0 

1.0 

1.2 

28.0 

4.0 

4.0 

10.6 

0. 

3.0 

2.6 

29.0 

3.0 

2.7 

21.9 

2.0 

1.0 

1.3 

28.7 

3.3 

3.3 

10.8 

2.7 

3.0 

2.7 

25.7 

6.3 

5.9 

8.7 

0. 

3.0 

2.9 

29.5 

2.5 

2.5 

13.5 

5.9 

2.0 

2.2 

26.5 

5.5 

5.5 

9.3 

0. 

3.0 

2.8 

29.0 

3.0 

3.0 

9.8 

5.5 

3.0 

3.0 

23.0 

4.0 

2.8 

20.  1 

0. 

i.O 

1.4 

25.8 

6.2 

6.2 

6.0 

2.8 

5.0 

4.3 

31.0 

1.0 

1.0 

18.7 

6.2 

1.0 

1.7 

31.0 

1.0 

1.0 

93.6 

0. 

0. 

0.3 

29.0 

3.0 

3.0 

17.4 

0. 

2.0 

1.7 

30.0 

2.0 

2.0 

26.5 

2.0 

1.0 

1.1 

31.0 

1.0 

1.0 

30.2 

2.0 

1.0 

1.0 

31.0 

1.0 

1.0 

30.8 

1.0 

1.0 

1.0 

31.0 

1.0 

1.0 

31.0 

1.0 

1.0 

1.0 

30.0 

2.0 

2.0 

30.0 

0. 

1.0 

1.0 

99 


QS 

X 

P 

Y 

Z 

CP 

30.  0 

2.0 

2.0 

30.0 

1.0 

1.0 

1.0 

23.0 

4.0 

3.6 

15.6 

0. 

2.0 

1.3 

23.6 

3.4 

2.3 

24.7 

1.6 

1.0 

1.2 

27.9 

4.  1 

4.  1 

10.  6 

2.3 

3.0 

2.6 

29.0 

3.0 

1.1 

55.  1 

1.1 

0. 

0.5 

23.  1 

3.9 

3.4 

16.5 

1.1 

2.0 

1.7 

30.5 

1.5 

1.5 

26.  7 

3.4 

1.0 

1.1 

28.5 

3.5 

3.5 

15.6 

0. 

2.0 

1.8 

29.0 

3.0 

2.  3 

24.9 

1.5 

1.0 

1.2 

27.3 

4.7 

3.  7 

14.9 

0.3 

2.0 

1.8 

23.0 

7.0 

5.5 

9.0 

0.7 

3.0 

2.3 

26.  5 

5.5 

4.3 

12.3 

3.5 

2.0 

2.2 
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I 


OR 

QS 


QSs.2. 

QR=30. 

CPs6. 

RLPHR=.8 

QLRH0Rs2. 


SIMULATION 


RESULTS  ***** 


INITIAL  CONDITIONS 

QS-50.^  QR-  5.0  CP-  5.0  ALPUA-O.UO  LAUDA-  5.0 


QS 

QR 

X 

P 

Y 

Z 

CP 

50.0 

5.0 

5.0 

10.0 

0. 

0. 

5.0 

53.0 

2.0 

2.0 

20.4 

5.0 

2.0 

2.6 

44.0 

11.0 

3. 1 

5.7 

0. 

9.0 

7.7 

45.  1 

9.9 

9.9 

3.  1 

6.  1 

5.0 

5.5 

50.0 

5.0 

5.0 

13.5 

6.9 

2.0 

2.7 

46.0 

9.0 

7.8 

5.9 

5.0 

9.0 

7.7 

50.3 

4.2 

4.2 

12.9 

7.8 

3.0 

3.9 

50.0 

5.0 

5.0 

10.4 

4.2 

5.0 

4.8 

44.0 

11.0 

11.0 

4.5 

5.0 

11.0 

9.8 

44.0 

11.0 

9.4 

5.3 

8.0 

3.0 

8.4 

51.4 

3.6 

3.6 

20.8 

8.4 

1.0 

2.5 

47.0 

3.0 

8.0 

7.7 

2.6 

7.0 

6.  1 

51.0 

4.0 

4.0 

11.5 

8.0 

4.0 

4.4 

47.0 

8.0 

8.0 

8.3 

2.0 

6.0 

5.7 

49.0 

6.0 

6.0 

11.3 

6.0 

4.0 

4.3 

48.0 

7.0 

7.0 

8.3 

5.0 

6.0 

5.7 

49.0 

6.0 

6.0 

9.5 

6.0 

5.0 

5.  1 

49.0 

6.0 

6.0 

14.3 

3.0 

3.0 

3.4 

51.0 

4.0 
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3.0 

18.9 

0. 

3.0 

2.8 

50.0 
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5.0 

3.9 

49.0 
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4.0 

2.4 

42.7 

0. 

1.0 

1.2 

49.4 

5.6 

5.2 

19.0 

2.4 
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0.7 

46.4 

8.6 

6.7 

13.8 

0.4 

6.0 

3.4 

52. 1 

2.9 

2.9 

24.0 

6.7 

1.0 

2.2 

50. 1 

4.9 

4.2 

24.0 

0. 

2.0 

2.  1 

50.3 

4.7 

4.7 

16.5 

4.2 

4.0 

3.0 

54.0 

1.0 

1.0 

26.7 

4.7 

1.0 

2.0 

146 


'iS  c^R  X P Y z cr 

54.0  1.0  l.O  35.7  1.0  1.0  1.5 

50.0  5.0  5.0  lo.l  0.  4.0  2.8 

53.0  2.0  2.0  28.2  4.0  1.0  1.9 

51.0  4.0  4.0  17.4  2.0  4.0  2.9 

55.0  0.  0.  37.4  4.0  0.  1.5 

53.0  2.0  2.0  30.6  0.  2.0  1.7 

52.0  3.0  2.7  38.0  0.  1.0  1.4 

53.  7 1.3  1.3  45.4  2.7  1.0  1.2 

53.  7 1.3  1.2  90.8  0.  0.  0.6 

52.7  2.3  1.6  66.3  0.  1.0  O.o 

50.3  4.7  4.7  21.0  1.6  4.0  2.4 

48.0  7.0  5.4  17.8  0.7.  3.0  2.7 

50.4  4.6  4.6  21.5  4.4  2.0  2.3 

■ 51.0  4.0  2.3  43.4  0.6  0.  1.2 
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